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ABSTRACT 


Experiments are described in which the mean and unsteady 
loads were measured on a single blade of a model of the 
controllable-pitch propeller on the DD-963 Class Destroyer. 
The experiments were conducted behind a model of the DD-963 
hull under steady ahead operation, hull pitching motions, 
and simulated acceleration maneuvers. The experimental tech- 
niques are outlined and the dynamometer and data analysis sys- 
tem described. 

The results show that all significant loads except radial 
force are predominantly of hydrodynamic origin. The circum- 
ferential variation of all measured components of blade load- 
ing is primarily a once-per-revolution variation, with the 
variation following approximately the variation of the tan- 
gential wake velocity. 

For sinusoidal pitching of the hull with maximum pitch 
angle of 1.85 degrees and a simulated full scale frequency 
of 0.16 hertz, the peak-to-peak circumferential variation of 
measured forces and moments increased by approximately 50 per- 
cent over the values without hull pitching. 

For simulated operation during an acceleration maneuver, 
the circumferential variation of measured forces and moments 
varied approximately as the product of ship speed and propeller 
rotational speed. At no time during the simulated acceleration 
maneuvers were the circumferential variations of loads as large 
as during full power steady ahead operation. 

For steady ahead operation, circumferential variation of 
loading determined from the model experiments agreed fairly 
well with full-scale data, but was substantially larger than 
the theoretically calculated values. 

For all conditions evaluated, the results follow close to 
previously reported results of similar experiments on a model 
of the FF-1088. 


ADMINISTRATIVE INFORMATION 

The work reported herein was funded by the Naval Sea Systems Command 
(NAVSEA 033), Task Area S0379-SLOO1, Task 19977. The work was performed 
under David W. Taylor Naval Ship Research and Development Center (DINSRDC) 
Work Unit No. 1544-296. 

The English system of units was used in the original calculations 
presented in this report. Therefore, all data are presented in the 
English units. However, the International System (SI) of metric units are 


shown in the text in parentheses following the English units. 


INTRODUCTION 

Major naval ships powered with marine gas turbines and using 
controllable-pitch (CP) propellers for thrust reversal are currently being 
added to the Fleet. Ships with CP propellers include the DD-963 Class, 
the FFG-7 Class, and the DDG-47 Class. 

Accordingly, the Navy has been conducting a research and development 
(R&D) program to establish the technology for producing reliable CP pro- 
pellers with delivered power in the range of 35,000 to 40,000 horsepower 
(26,000 to 30,000 Wee As part of this program, CP propellers were 
installed on the U.S.S. PATTERSON (FF-1061) and U.S.S. BARBEY (FF-1088) 
with delivered power of 35,000 horsepower (26,100 kW). These installations 
were intended to demonstrate that CP propellers in this range of power had 
adequate reliability for application to ships with gas turbine prime 
movers. 

Because of the structural failure of the crank rings to which the 
blades of the CP propeller on the FF-1088 were bolted, R&D efforts were 
intensified. The program undertaken at DINSRDC included: 

1. Blade Loading of CP Propellers 

a. Model measurement and theoretical prediction of blade load- 
ing on CP propellers. 

b. Model and full-scale wake measurements and theoretical pre- 
dictions of wake. 

c. Full-scale measurements of forces, pressures, and strains in 
CP propeller components. 

2. Structural Design of CP Propeller Blade Attachments. 

3. Development of Materials for CP Propeller Systems. 

The current report presents the results of work conducted under Sec- 


tion la of the CP Propeller Research and Development Program, i.e., model 


Theses silie J.J. et al, "U.S. Navy Controllable Pitch Propeller Programs," 
presented at a Joint Session of the Chesapeake Section of the Society of 
Naval Architects and Marine Engineers and the Flagship Section of the 
American Society of Naval Engineers, Bethesda, Maryland (19 April 1977). 


measurement and theoretical prediction of blade loading of CP propellers. 

Work under the other sections of this program will be reported separately. 
The present report presents experimental results obtained on a model 

of the CP propeller on the DD-963 Class Destroyer. The results of similar 


experiments on a model of the FF-1088 were reported in Reference 2. 


BACKGROUND 

Extreme care must be taken to design the blades and pitch-changing 
mechanisms of high power CP propellers so that they possess adequate 
strength including consideration of yield and fatigue stresses. This 
requires an accurate estimate of the maximum time-average and alternating 
loads under all operating conditions. High time-average and alternating 
loads occur at steady full-power ahead conditions and during high-speed 
maneuvers including full-power crash astern, full-power crash ahead, and 
full-power turns. In addition, the influence of the seaway may substan- 
tially increase the time-average and alternating loads. At present there 
appears to be no confirmed technique whereby the pertinent loads can be 
predicted to the desired accuracy. Schwanecke and Weeeldeneae reviewed the 
factors affecting blade loading for propellers in general, and Rurseesiieey 
and Hawdon et ee discussed some of the factors peculiar to blade loading 


of CP propellers. 


Sogetaiil, R.J. et al, “Experimental Determination of Mean and Unsteady 
Loads on a Model CP Propeller Blade for Various Simulated Modes of Ship 
Operation," The Eleventh Symposium on Naval Hydrodynamics sponsored Jointly 
by the Office of Naval Research and University College London, Mechanical 
Engineering Publications Limited, London and New York, pp 789-823, 832-834, 
(April 1976); also "Experimental Unsteady and Mean Loads on a CP Propeller 
Blade of the FF-1088 for Simulated Modes of Operation," David Taylor Naval 
Ship Research and Development Center Report 76-0125, October 1976. 


Scehmemcekes H. and R. Wereldsma, "Strength of Propellers Considering 
Steady and Unsteady Shaft and Blade Forces, Stationary and Nonstationary 
Environmental Conditions," Proceedings of the Thirteenth International 
Towing Tank Conference, Report of the Propeller Committee, Appendix 2b, 
Vol. 2, pp 495-526 (1972). 


sRucetelinn A.A., "Hydrodynamics of Controllable Pitch Propellers," 
Shipbuilding Publishing House, Leningrad (1968). 


5 

Hawdon, L. et al, "The Analysis of Controllable-Pitch Propeller Char- 
acteristics at Off-Design Conditions," Transactions of the Institute of 
Marine Engineers, Vol. 88, Series A, Part 4, pp 162-184 (1976). 


3 


Near the self-propulsion point in calm water, the time-average loads 
can probably be calculated with reasonable accuracy. However, even at 
these conditions, the variation of loads with blade angular position ap- 
parently cannot be calculated with high accuracy. Various techniques, 
including quasi-steady procedures, stripwise unstead procedures, and 
methods based on unsteady lifting surface theory, have been proposed for 
calculating the unsteady loading arising from the circumferential varia- 
tion in the inflow <gallaeiiey 5° However, all of these procedures require 
knowledge of the flow patterns (wake profile) in the propeller disk. In 
current practice, the wake profile is measured in the plane of the propel- 
ler behind the model hull with the propeller removed. For high-speed dis- 
placement ships of the type under consideration in this report, these 


results are usually extrapolated to full scale without making allowance for 


Srercusien, J.E., "The Development of Numerical Methods for the Computa- 
tion of Unsteady Propeller Forces," Presented to the Symposium on Hydro- 
dynamics of Ship and Off-Shore Propulsion Systems, Oslo, Norway (March 
UD TD) 


Tteon T. et al, "Calculation of Unsteady Propeller Forces by Lifting 
Surface Theory," Presented to the Symposium on Hydrodynamics of Ship and 
Off-Shore Propulsion Systems, Oslo, Norway (March 1977). 


eroddy: R.F., “A New Method for the Calculation of Unsteady Forces on 
a Marine Propeller," Presented to the Chesapeake Section of the Society 
of Naval Architects and Marine Engineers, Washington, D.C. (February 1977). 


aon Gent, W., "Unsteady Lifting Surface Theory for Ships Screws: Deri- 
vation and Numerical Treatment of Integral Equations," Journal of Ship 
Research, Vol. 19, No. 4, pp 243-253 (December 1975). 


psehwaneekes H., "Comparative Calculations on Unsteady Propeller Blade 
Forces,'' Proceedings of the Fourteenth International Towing Tank Confer- 
ence, Report of Propeller Committee, Appendix A, Vol. 3, pp 357-397 (1975). 


teeeciltim, J.P., "Propeller Excitation Theory," Proceedings of the 
Thirteenth International Towing Tank Conference, Report of the Propeller 
Committee, Appendix 2c, Vol. 2, pp 527-540 (1972). 


a rerreiele R.J. and M.L. Miller, "Unsteady Propeller Loading - Measure- 
ment, Correlation with Theory, and Parametric Study," Naval Ship Research 
and Development Center Report 2625 (October 1968). 


(1) the change in Reynolds number and the corresponding reduction in 
relative boundary layer thickness and (2) the effect of the propeller suc- 
tion on the boundary layer and thereby the wake pattern in propeller disk. 
Although these two factors may be important for full-form ships such as 
cargo ships, they are probably not important for high-speed transom stern 
ships of the type under consideration in this report. 

Existing measurements which give information on unsteady blade load- 
ing include: 

1. Measurements of strain on the blades of the model propellers or 
full-scale propellers. However, some calculations and assumptions are 
required to convert measured strains into loads. Published data of this 
type have been summarized by Moyne 

2. Measurement of bearing (shaft) forces and moments on model pro- 
pellers operating in wakes generated by model hulls or wire grid screens. 
However, these experiments give information on only some components of 
blade loading and on only those harmonics of shaft rotational speed cor- 
responding to nZ-1, nZ, and nZ+l, where n is an integer and Z is the 
number of blades. Measurements of this nature have been conducted by many 
investigators, as summarized by Baseline and Weceldenatue 

3. Measurements of forces and moments on individual blades of model 
propellers operating in wakes generated by model hulls or wire grid 
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screens. Measurements behind model hulls have been made by Huse and 
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Micrel measurements behind screens have been made by Hawdon et al, 
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and Rees caution measurements in inclined flow have been made by Albrecht and 
Suinnlies Predinenenilleoe and Raccundl. “ere measurements on partially sub- 
merged propellers have been made by DOBeveas 

Experiments in wakes generated by screens are advantageous for 
evaluating the ability of a procedure to calculate the loading for a given 
wake since for this case, the propeller apparently does not influence the 
wake pattern. Although some good correlation has apparently been obtained 
between analytical predictions and unsteady bearing forces measured be- 


Tak aL 


hind wire grid screens, correlation has been rather inconsistent 


between analytically predicted unsteady blade loads, or resulting 


; S aL 21 
strains, and measured blade loads, or strains. 20 


The mechanism by which the seaway influences the mean and unsteady 
blade loads is complex. Factors include the increased mean propeller 
loading due to increased hull resistance and the increased unsteady load- 
ing resulting from the influence of the free surface and modification of 
the flow pattern into the propeller disk. This flow pattern is influenced 
by (1) direct trochoidal velocities from the ocean waves, (2) relative 
velocities of the propeller due to ship motions, and (3) modification of 
the hull wake pattern due to the seaway and ship motions. Procedures for 


calculating the loads in a seaway are much less refined than for steady 
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operation in calm water. mayeatelag gives a good review of the mechanisms 
and procedures for predicting the effect of the seaway on bearing forces 
which, in principle, also apply to unsteady loading on an individual 
blade. Keil et ata and Watanabe et a a present strain measurements 
on the blades of full-scale propellers in both calm and rough seas. 

Apparently no rational analytical procedures are available for 
accurately calculating the time-average loads per revolution or the un- 
steady loads including variation with blade angular position during crash- 
ahead or crash-astern maneuvers. These loads may depend on many factors 
including the time rate of change of propeller pitch p (for CP propel- 
lers), time rate of change of rotational speed n, time rate of change of 
ship speed V, propeller blade-section stall, cavitation, ventilation, 
flow separation from the hull, and large interactions between the propel- 
ler and the hull. Some of these factors are discussed and considered by 
Hawdon et als? 

For turns, the factors affecting the time-average loads per revolu- 
tion and the unsteady loads are somewhat the same as those affecting the 
loads under crash-ahead and crash-astern conditions except that for turns, 
there is a relatively large drift angle of the flow into the propeller. 
This drift angle tends to increase the circumferential nonuniformity of 
the flow into the propeller, thereby increasing the unsteady loading. 
However, this circumferential nonuniformity of the inflow tends to be off- 
set by the lower values of ship speed and propeller rotational speed in 


turns compared to steady ahead operation. 


Prior to the present R&D investigation, no experimental measurements 


existed to the authors' knowledge which showed the time-average loads and 
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circumferential variation of loads with blade angular position on CP pro- 
pellers behind a hull under a wide range of operating conditions. An ex- 
perimental program was therefore undertaken to measure the six components 
of loading (Figure 1)* on model CP propellers operating behind model hulls. 
The initial experiments were conducted on a single-screw ship, namely, a 
model of the FF-1088. These results were reported winreriomoliy The sec- 
ond set of experiments were conducted on a twin-screw ship, namely, a 
model of the DD-963 Class. These results are presented in the current 
report. 

For the DD-963 Class, the experimental conditions included 
(1) steady-ahead operation near the self-propulsion point, (2) steady- 
ahead operation near the self-propulsion point with forced dynamic pitch- 
ing of the model hull, and (3) simulated acceleration operation. 

Results for the steady ahead operation were correlated with predic- 
tions based on unsteady lifting surface theory as developed by Tsakonas 
et ail”? and with the quasi-steady method of WeCercttinnn@° and with strains 
measured on the full-scale propeller. 

Blade loading measurements were made on the propeller on the star- 
board shaft since this shaft had a larger rake angle than the port shaft. 
The propellers used in these experiments were DINSRDC propellers 4660 
(right hand rotation on port shaft) and 4661 (left hand rotation on star- 
board shaft), which were made of aluminum; see Figure 2 and Table 1.** 
The hull of the DD-963 Class was represented by DINSRDC model hull 
5265-1B; see Figure 3. 


*Figures are presented following the section on acknowledgments. 


**kThe tables are presented following the figures. 
2 ecikones, S. et al, "An Exact Linear Lifting Surface Theory for 
Marine Propeller in a Nonuniform Flow Field," Journal of Ship Research, 
Vol. 17, No. 4, pp 196-207 (December 1974). 

6 
AO eve tee. J.H., "On the Calculation of Thrust and Torque Fluctuations 
of Propellers in Nonuniform Wake Flow," David Taylor Model Basin Report 
1533 (October 1961). 


EXPERIMENTAL TECHNIQUE 


FACILITY AND DYNAMOMETRY 

All experiments were conducted on DINSRDC Carriage I using basically 
the same dynamometry and hardware as previously described in Reference 2. 

The port propeller, on which blade loads were not measured was driv-— 
en from inside the model hull as would be the case in a self propulsion 
experiment. The propeller rotational speed, which could be controlled 
independently of the starboard propeller, was measured via a toothed- 
pickup and recorded on a digital voltmeter. The time-averaged thrust and 
torque were measured for selected runs by a transmission dynamometer. 

The starboard propeller, on which blade loads were measured, was 
located in its proper position relative to the model hull but was isolated 
from the hull and driven from downstream (see Figure 4). This downstream 
drive system was necessary in order to obtain the required characteristics 
of the system for measuring unsteady loading. The general criteria for 
the design of an unsteady force measuring system are: : 

1. The support structure of the force measuring system should be 
soft mounted and possess a large mass to eliminate transmission of extra- 
neous vibration to the system. 

2. The natural frequency of the system should be well above the 
highest frequency of the quantities to be measured (to avoid phase shift 
and amplification of the signal). 

3. The system response in the force magnitude range should be suf- 
ficiently large to be measurable (sensitivity). 

4. The system should be free of interaction, that is, each measur- 
ing element should respond only to that force or moment which it is in- 
tended to measure. 

These four major aims are not complementary. The high natural fre- 
quency requires a stiff, rigid system whereas high sensitivity requires 
an elastic, soft system. The necessary compromise results in some inter- 
action between the force-measuring elements. 

Criterion 1 dictated that a massive flywheel be used, and Criterion 
2 dictated that this flywheel be connected: to the sensing elements (locat-— 


ed inside the propeller hub) by a short thick shaft. Therefore, because 


of the geometry of the hull and shafting of the configuration under eval- 
uation, it was not feasible to achieve both these criteria with an up- 
stream drive system from inside the model hull. Criteria 1 and 2 con- 
trolled the minimum allowable beam and draft of the downstream body and 
the maximum allowable clearance from the bow of the downstream body to the 
propeller. Although the downstream body may exert some influence on the 
flow into the propeller, that location was considered necessary in order 
to meet these measuring criteria. The influence of the downstream body 
on the flow into the propeller is discussed in the section on experimental 
results. 

The drive and mounting system was basically the same as that used 
in the DINSRDC BASS dynamometer which has been described by Beaadedeaa 
Utilized from this dynamometer were the propeller (tail) shaft, drive 
shaft with flywheel, belt-type (quiet) transmission, and sliprings. 
Power to rotate the propeller was supplied by a d-c permanent-magnet 
servomotor capable of delivering up to 33 foot-pounds (45 N-m) of torque. 
The electrical power to this motor was delivered through a precision 
solid-state motor controller so that the shaft revolution rate could be 
controlled very accurately and held over the wide range of propeller 
torque loadings required for some of the experimental conditions. Mounted 
on the propeller shaft was a digital encoder that generated electrical 
pulses as a function of shaft angular position. Two types of pulses were 
generated: a single pulse per revolution and a multipulse per revolution 
(90 equally spaced pulses for the current experiment). The single pulse 
was syncronized with the reference line of the instrumented propeller 
blade. The pulses generated by this encoder are accurate to within 0.01 
degree. 

The downstream body which housed the drive system was basically that 
used by eibag but modified to allow deeper submergence and an inclined 
shaft angle. Both the body housing (the drive system was soft mounted to 


this body) and the model hull were rigidly attached to a pitch-heave 


or less, J.H., "Static and Dynamic Calibration of Propeller Model 


Fluctuating Force Balances," David Taylor Model Basin Report 2350 (March 
1967); also Technologia Naval, Vol. 1, pp 48-74 (January 1968). 
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oscillator which, in turn, was rigidly mounted on the towing carriage. 
This arrangement enabled the model hull and the drive system to be dy- 
namically pitched together while maintaining independent support from one 
another. 

The sensing elements were flexures to which were bonded high- 
sensitivity, semiconductor strain-gage bridges. The design of these 
flexures has been described by Dobayeae There were three flexures, each 
of which measured two components of blade loading. Flexure 1 measured 
components Ee and Hho Flexure 2 measured components Be and M> and Flexure 
3 measured components Ee and M, (Figures 1, 5, and 6). An arrangement of 
three separate flexures rather than one to measure all components of blade 
loading was adopted Because it appeared to result in higher natural fre- 
quencies (Criterion 1), higher sensitivities (Criterion 3), and lower 
interactions (Criterion 4), than would have resulted had a single flexure 
been used. 

The flexures were mounted inside a propeller hub which was specifi- 
cally designed for these experiments (Figure 6). Only one flexure could 
be mounted at a time, because of space limitations, and this necessitated 
duplicate runs, as discussed later in the section on experimental condi- 
tions and procedures. 

The strain-gage bridges were excited by a common d.c. voltage 
source, transmitted through the sliprings on the propeller shaft. The 
constant-current excitation used by DObayae was not employed in the pre- 
sent experiment because it appeared to be too sensitive to temperature. 

The voltage output from the flexures (due to blade loading) was 
transmitted through the sliprings to individual amplifiers (NEFF 119-121). 
These amplifiers utilized field effect transistors to produce an extremely 
high input-impedance (100 megohms, minimum). This high impedance essen- 
tially eliminated slipring noise to the amplifier. The voltage signals 
were transferred across the sliprings in the presence of only a small 
amount of noise-producing current. The amplifiers used here had zero- 
phase shift qualities in the d.c. to 20 kilohertz range. They were chopper- 
stabilized to enable both the steady and unsteady signals to be recorded 
simultaneously. This signal-conditioning system was essentially the same 


as that used by Dope. 


iil 


The signals were then digitized and analyzed by using a Model 70 
Interdata Digital Computer, and then stored in digital form on a nine- 
track magnetic tape. The on-line analysis of the data is discussed in the 


section on data acquisition and analysis. 


CALIBRATION 

Prior to the experiment, each flexure was statically calibrated in 
air to establish flexure sensitivities, interactions, and linearity over 
the loading range of interest. These calibrations were conducted with the 
flexures mounted in a calibration stand, with the flexure electrical ca- 
bles connected through the flywheel and drive assembly as in the experi- 
ment. Each flexure was subjected to independently controlled forces in 
the axial, transverse, and radial directions (i.e., Flo ee and Fo re- 
spectively) and to independently controlled moments about the axial, 
transverse, and radial directions (i.e., M, Bye and M,> respectively); 
see Figure l. 

The static calibration showed that all flexures had a linear re- 
sponse over the load range of interest. Table 2 shows the interaction 
matrix. These calibrations indicated that all flexures had good sensi- 
tivity except Be whose sensitivity was slightly lower than desirable. 

The interactions were small except for the effect of M, on ae The infe- 
rior characteristics of the ER flexure is not considered a serious short- 
coming since Ek arises primarily from centrifugal loading and can be 
analytically calculated. In addition, no significant variation of ES 
with blade angular position was anticipated. Flexure 3, which measured 
F. and M,> was further evaluated by correlation of air-spin experiments 
with analytically calculated results, as discussed later. The interac-— 
tions were taken into consideration during data analysis. 

The flexures used in this experiment had been dynamically cali- 
brated by easy” to determine the frequency range over which unsteady 
forces and moments could be reliably measured. In this procedure, an 
electromagnetic shaker in air was used to apply a relatively constant, 
maximum amplitude, variable-frequency force or moment-excitation in all 
six-component directions to all six flexure elements. The force or mo- 


ment amplitude imposed by the shaker was monitored through an extremely 
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light-weight, strain-gaged single flexure element. The measured lowest. 


natural frequencies of the three flexures in air were as follows: 


Frequency (hertz) Mode 


Flexure 1 550 M, 
Flexure 2 450 MO 
Flexure 3 282 M, 


The measured amplification factor (ratio of output amplitude to in- 


put amplitude) and phase shift for all three flexures was as follows: 


Frequency Range (hertz) 0 to 60 60 to 120 
Phase Shift (degrees) 0 to 0.05 0.05 to 0.15 
Amplification Factor 1.00 1.00 to 1.05 


In the previous eeperimentd the natural frequency of the flexures 
was found to be reduced by 55 percent when submerged with blades attached. 
As a result, it was concluded that the flexures had a "true" dynamic re- 
sponse up to at least the third harmonic and no greater than a five per- 
cent amplification up to the sixth harmonic. Because the blades used on 
the present experiment were lighter and smaller than those used on the 
previous experiment, it was assumed that the natural frequency of the 
flexures would be reduced to a lesser extent when submerged with blades, 
so that the dynamic amplification would be less. This assumption was 
supported by the increase in frequency of the extraneous signals appearing 
in the unfiltered experimental data as discussed in the experimental re- 
sults section. 

The propeller shaft drive and soft-mount support system were dynam- 
namically loaded in the vertical, longitudinal, and transverse directions 
to obtain the lowest natural frequencies of the system. The natural fre- 


quencies of the system in air were found to be: 


Mode Natural Frequency (hertz) 
Vertical bending W252 
Horizontal bending 6.0 
Axial 4.6 


The support system had a low resonant range; however, the soft- 


mount system was specifically designed to prevent towing-carriage 
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oscillation (with the resonance at 100 to 200 hertz) from being trans- 
mitted to the blade flexures. Based on the measured resonance, it is 
concluded that the soft-mount system should successfully meet this ob- 
jective. Although some resonances were close to the propeller rotational 
speed for some experimental conditions, it was considered more desirable 


to isolate the system from towing-carriage vibration. 


EXPERIMENTAL CONDITIONS AND PROCEDURES 

Experiments were conducted at several conditions including steady- 
ahead operation, simulated pitching of the hull, and simulated accelera- 
tion. All conditions were run with the model hull rigidly attached to its 
support, with no freedom to sink or trim, and with essentially equal rota- 
tional speed on the port and starboard propellers. 

The steady-ahead condition is defined in Tables 3 and 4. The simu- 
lated full scale ship speed and propeller rotational speed for this condi- 
tion were determined from model self-propulsion data* at simulated dis- 
placement of 7,800 tons (7,920 tonnes) including corrections for wind drag 
at zero true wind and a three-percent margin on effective power with 
Cy = 0.0005. 

The trim and draft at this speed were obtained from Reference 28. 
These had been determined by setting the specified still water trim (even 
keel) and draft (19.5 feet (6.40 m) full-scale equivalent), attaching 
the model to the carriage so that it was free to trim and sink, running at 
the specified speed, and locking the model at this equilibrium trim and 
draft. The equilibrium sinkage was 0.5 feet (15.4 cm) at the bow and 3.0 
feet (98.4 cm) at the stern. 

Runs simulating hull pitching were conducted at the same conditions 
as the steady-ahead run, except that the hull pitch was varied. Two types 
of runs were conducted: (1) quasi-steady simulation in which the hull 


pitch angle ~ was set at various fixed positions, and (2) unsteady 


*DTNSRDC experiments 21 and 22 on Model 5265-1B. 

eve W.G., "The Effect of Speed on the Wake in Way of the Propeller 
Plane for the DD-963 Class Destroyer Represented by Model 5265-1B," David 
Taylor Naval Ship Research and Development Center Report SPD-311-37 (July 
UGS) ¢ 
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simulation in which was varied sinusoidally with time. For the quasi- 
steady simulation, runs were conducted at five different values of yw, from 
1.85 degrees bow up from the calm water equilibrium VW=V ou) to 1.85 de- 
grees, bow down from vow (Tables 3 and 4). For the unsteady pitch simu- 
lation, the value of ~ was varied sinusoidally about vow with an ampli- 
tude of 1.85 degrees and a frequency of 0.8 hertz.* The selected scaled 
amplitude and frequency were within the predicted response characteristics 
of the DD-963. All runs were conducted in calm water; therefore, the re- 
sponse of the hull to the seaway was simulated but the seaway was not 
simulated. 

Acceleration runs were conducted based on analytical dynamic simu- 
lation studies of the DD-963.-> The experimental conditions followed run 
7501061 of Reference 29, which was an acceleration from 8.7 knots to full 
power (see Tables 3 and 4 and Figure 7). Trim and displacement were fixed 
at the values corresponding to the self-propulsion condition (Condition 1 
of Table 3). Two types of runs were conducted: (1) quasi-steady runs in 
which all quantities including model speed V, rotational speed n, and pro- 
peller pitch P were held constant (V=n=P=0), and (2) unsteady runs in 
which V was varied with time but n and P were held constant (V>0, n=P=0). 
For the quasi-steady simulation, runs were conducted at five different 
combinations of V, n, and P. The conditions for each run represent the 
conditions at one instant of time during a "true'' acceleration in which 
V, n, and P vary with time. Thus, one "true" acceleration run is repre- 
sented by five steady runs which do not simulate the time rate of change 
of V, n, and P. For the unsteady simulation, runs were conducted at the 
same five combinations of fixed n and P as used for the quasi-steady simu- 
lation, and V was varied with time (the same variation was used for each 
run) representing an acceleration of the model hull (Figure 7). For each 
of these runs, data are of interest only near that value of V which oc- 


curred concurrently with the fixed values of n and P during the "true" 


x 
Full scale equivalent frequency is 0.16 hertz. 
Cea C.J. and T.R. Harper, "Propulsion Dynamics Simulation of the 


DD-963 Class Destroyer," Propulsion Dynamics, Inc., Report 74R1B (January 
OID) he 
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acceleration (V>0, n#0, P#0). Thus, one "true" acceleration run is rep- 
resented by five runs which simulate the proper time rate of change of V 
but not the proper time rate of change of n and P. The quasi-steady and 
unsteady acceleration simulations were for the same conditions, the only 
difference being that V=0 for the quasi-steady simulation whereas V>0 for 
the unsteady simulation. In general, P varies with time during a "true" 
acceleration run; however, for the acceleration run under simulation here, 
P was constant throughout the portion of the run simulated. 

For the unsteady acceleration runs, the carriage speed was manually 
varied with time in a carefully controlled manner. This was achieved with 
the aid of an inked pen on a two-dimensional Cartesian plotter. In one 
direction, the pen was controlled so that it moved linearly with time, 
and in the orthogonal direction, it was controlled so that it varied with 
the instantaneous carriage speed. When an acceleration maneuver was to be 
executed, the switch moving the pen with time was turned on and the car- 
riage operator manually varied the carriage speed so that the inked pen 
followed a prescribed velocity versus time curve. 

As discussed earlier, each of the three load-sensing flexures 
measured only two components of blade loading. Therefore, each of the 
experimental conditions described in Table 3 was run with each of the 
three blade loading flexures. 

The blade pitch was set by using a Sheffield Cordax 300 measuring 
machine. In order to change either the blade pitch or the flexure, the 
propeller had to be removed from the drive system. 

Air-spin experiments were conducted with all three flexures over a 
range of rotational speeds in order to isolate the effects of centrifugal 
and gravitational loading from hydrodynamic loading. Supplemental experi- 
ments were conducted to assess the influence of the downstream dynamo- 
meter boat on the flow in the propeller plane. These supplemental exper- 
iments consisted of wake surveys in the propeller plane at the self- 
propulsion point (Condition 1 in Table 3) with and without the downstream 
body, but without the propeller. These wake surveys yielded a direct 


measure of the change in the velocity distribution through the propeller 
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disk attributable to the downstream body. The details of these wake sur- 


veys will be reported in a future DINSRDC report.* 


DATA ACQUISITION AND ANALYSIS 

Data were collected, stored, and analyzed on-line by using a Model 
70 Interdata Digital Computer. A special-purpose computer program was 
written with options for analyzing each of the three basic types of runs: 
(1) steady ahead, (2) dynamic hull pitching, and (3) unsteady acceleration. 
These types of runs have already been discussed in detail. 

The program allowed the propeller blade force and moment data to be 
sampled and stored on magnetic tape as a function of shaft position. 
Sampling was triggered by external pulses generated by a Baldwin Digital 
Encoder mounted on the propeller shaft, as discussed earlier. Pulses were 
generated as a function of shaft angular position; hence, the sampling of 
blade force and moment data was related to shaft position. There were two 
outputs from the shaft encoder; a single pulse per revolution and multi- 
ile (OO pulses per revolution for the current experiments). 

When the experimental condition was achieved, the computer operator 
initiated the data collection cycle. The program "waited" until the sin- 
gle pulse occurred; when the single pulse occurred, the computer again 
"waited" for the occurrence of the first following pulse of the 90 pulses; 
data were then sampled for all channels through an analog-to-digital con- 
verter and stored in computer memory. This process was repeated for 180 
pulses, or two shaft revolutions. At the same time, the program "read" 
two frequency counters into core memory which measured model velocity V 
and propeller rotational speed n. The values of V and n were measured by 
counting the pulses from geared wheels attached to the towing carriage 
drive system and to the propeller shaft, respectively. The values of V and 
n were averaged over two shaft revolutions. Thus, there was an average V 
and n corresponding to each pair of two consecutive revolutions. 

After two revolutions of data were sampled and stored in core memory, 
the data were transmitted from core to a nine-track digital tape recorder. 


The transfer time was small and no pulses were missed during the transfer. 


* 
These wake surveys were conducted by R.F. Roddy, DINSRDC Code 1524. 
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The data collection cycle proceeded continuously until the operator dis-— 
engaged the computer. The sampling procedure was the same for all types 
of experimental conditions, and at the completion of an experimental run, 
all data were stored on magnetic tape and were available for analysis im- 
mediately or at any later time. For the analysis, the computer operator 
selected the appropriate option of the program depending on the type of 
run, i.e., (1) steady ahead, (2) dynamic hull pitching, or (3) unsteady 
acceleration. 

The appropriate calibration factors were stored in the computer and 
considered in the analysis. However, since only two of the six components 
of blade loading were measured during a given run, the interactions between 
the various loading components could not be considered during the on-line 
analysis. The interactions were taken into account later after measure- 
ments were completed with all three flexures for a given condition. 

For the steady ahead condition, blade force and moment data at each 
4-degree increment of blade angular position were averaged over the number 
of cycles recorded (usually over more than 200 cycles). Spurious data not 
related to shaft position are averaged out by this method. An harmonic 
analysis was then performed on the average wave forms of the blade loading 
components. This gave the amplitude and phase of the first 16 harmonics. 

For the dynamic pitch runs, the hull pitch angle varied sinusoidally 
with a frequency of 0.8 hertz. A position potentiometer translated bow 
vertical displacement into hull pitch angle, and this was read into the 
computer in the same manner as blade loading components. During dynamic 
pitching, the shaft rotated independent of the pitch oscillator. During 
a single propeller revolution, 90 pitch positions were measured. Thus, 
to correlate pitch angle position and revolution, an average pitch must be 
taken over each revolution. 

Fourteen dynamic pitch angle positions were selected for analysis. 
These were characterized by pitch angle ) and the sign of the time rate of 
change of pitch angle j. The computer calculated the average ~ and sign 
of corresponding to each propeller revolution. Based on these calcu- 
lated average values of w and sign of }, each propeller revolution was 
either placed in a suitable hull-pitch angle category or discarded if its 


average ) fell outside the tolerance band of all the 14 specified values 
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of ~. Several passes down the towing tank were required in order to obtain 
a sufficient number of samples. After all the data had been sorted based 
on ~, sign of }, and tolerance, the cycles for each combination of jp, and 
sign of ~ were analyzed in exactly the same manner as the data for the 
steady-ahead condition at fixed jp. 

For unsteady-acceleration runs, the model speed V varied with time 
t. During an acceleration run, data, including a measure of V, were sam- 
pled and stored in the same manner as for the steady-ahead runs. 

Five values of V were specified for analysis. For each acceleration 
run and for each specified V, the computer selected the propeller revolu- 
tion which had the average value of measured V nearest to the specified V. 
However, because only one revolution at each specified velocity was ob- 
tained for a single acceleration run, each such run was repeated five 
times. This yielded five revolutions at each specified velocity. All the 
cycles for each specified V were then analyzed in exactly the same manner 
as the data for the steady-ahead conditions. 

Thus, the on-line analysis system yielded average wave forms and 
harmonic analyses of the average wave forms for steady-ahead conditions, 
for specified conditions of ~, sign of W during the dynamic pitch cycle, 
and for specified velocities V during the acceleration operation. However, 
these on-line results are preliminary because: 

1. They do not consider the interactions between the various load 
components. These interactions were determined during the static cali- 
bration of the flexures. 

2. They include the complete measured signals with no filtering. 

As discussed in the section on experimental results, some extraneous sig- 
nals near the natural frequency of the flexure being used appeared to be 
superimposed on the signals generated by blade loading. 

3. They include the effect of centrifugal and gravitational load- 
ing on the aluminum model propeller. The centrifugal and gravitational 
components of loading were measured separately during air-spin experiments, 
as discussed earlier. 

4. They do not have any corrections for the influence of the dyna- 


mometer boat. These corrections are discussed later. 
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5. The bending moments were calculated about the radial location of 
the strain gages on the flexures, rather than about the shaft axis or some 
desired radius on the blade. 

Final analyses were conducted after completion of the experiment to 
consider interactions, to filter out extraneous high frequency noise, to 
isolate hydrodynamic loading from centrifugal and gravitational loading, 
to correct for the dynamometer boat, and to resolve bending moments as 
desired. These analyses were conducted using a Control Data Corporation 
(CDC) 6700 Computer. 

For each condition, the average wave form for each of the six load- 
ing components was multiplied by the inverse of the calibration matrix 


given in Table 2. 
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This matrix multiplication was performed at 4-degree increments of blade 
angular position. 

An harmonic analysis was then performed on the signals corrected for 
the interactions. Based on an harmonic analysis of the wake in the pro- 
peller wile, it was judged that there should be no significant loading 
of hydrodynamic origin at frequencies above ten times shaft frequency. 
Therefore, the wave form was then reconstructed by using the first ten 
harmonics of shaft frequency. This reconstruction using only the first 
ten harmonics had the same effect as filtering out all frequencies above 


ten times shaft frequency. 


Cecmtincss D.E., "Numerical Prediction of Propeller Characteristics," 
Journal of Ship Research, Vol. 17, No. 1, pp 12+18 (March 1973) 
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Corrections were made to the mean values of the measured loading 
components to account for centrifugal loads and the influence of the dy- 
namometer boat, and to the first harmonic of the measured loading compo- 
nents to account for gravitational loads. The derivation of these correc- 
tions is discussed later. 

From the known values of the three measured force components and 
three measured moment components, the values of the bending moments about 
the shaft centerline and bending moments normal to the nose-tail line at 
the 0.3 and 0.4 radii were calculated. These bending moments were calcu- 
lated at every 4 degrees of blade angular position, and harmonically ana- 
lyzed. The wave form was reconstructed by using the first 10 harmonics of 
blade angular position, in exactly the same manner as was used for the 
other components of blade loading. 

Plots of the data were generated by the CDC computer system using a 


-Calcomp Plotter. 


ACCURACY 

During the experiments for steady operation, V=0, and dynamic 
pitching, W#0, where many revolutions of data were averaged during a sin- 
gle run, the standard deviations of speed V, rotational speed n, forces, 
and moments were computed, assuming the distribution in these variables at 
a given condition follows the normal probability distribution. For forces 
and moments, the standard deviation was calculated at every increment of 
blade angular position at which forces and moments were recorded. An 
error band around the data mean was then represented using the standard 
deviation multiplied by a factor dependent on the confidence level chosen. 
For the present analysis, the factor of 1.96 was selected which corre- 
sponds to a confidendence level of 95 percent. A confidence level of 95 
percent indicates a confidence (or probability) that 95 percent of the 
data considered falls within the error band. For a given run the average 
error (95 percent confidence band) in model speed V was approximately 
+0.005 foot per second (1.6 mm/s), while the error in rotational speed n 
was less than 0.001 revolution per second. The very low error in n re- 
sulted from the use of a precision solid-state motor controller as dis- 


cussed in the section on facility and dynamometry. 
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For a given steady run (vV=0, #0) the error (95 percent confidence 
band) in measured forces and moments, fluctuated from +5 to +10 percent of 
the mean signal, depending on the circumferential blade position. Figure 
8 demonstrates the variation in error in one revolution of the uncorrected, 
raw Bas signal at the full power condition. This figure represents the 
general trend of all the force and moment components measured. 

Besides the fluctuation in signals occurring in a given run, the 
overall accuracy of the data is dependent on the repeatability from one 
run to the next. An effort was made to set experimental conditions iden- 
tically on repeat runs; however, the propeller rotational speed and model 
velocity were set by hand, so some variation was unavoidable. Table 5 
demonstrates the variation in the measured experimental conditions and the 
raw data for the Ee component for 11 repeat runs. The variation in the 
mean force was +4 percent over all the runs, but on a given day the vari- 
ation averaged +2 percent. The same trend can be observed in rotational 
speed, model velocity and the harmonic force components. This day-to-day 
variation could be due to different operators setting the experimental 
conditions, slight variations in the draft of the model, and variations in 
the gain of the sensing electronics. The variations shown for a are 
typical of all the measured force and moment components. 

For all experimental conditions the rotational speed of the port and 
starboard propellers were intended to be equal. However, some exploratory 
runs were conducted to determine whether the mean or unsteady loads, which 
were measured on the starboard propeller, were influenced by the rotational 
speed of the port propeller. At a fixed value of rotational speed on the 
starboard propeller ne» the rotational speed on the port propeller n_ was 
varied. The results showed that there was no measurable effect of n_ on 
the mean or unsteady loads in the region 0.95 ne £8, < 1.05 no: For all 
runs for which data are presented, 0.99 n. 5 < 1.01 ns; therefore, the 
results presented are not measurably influenced by inaccuracies in aye 

For the unsteady acceleration (V>0), the average of the five values 
of V and n for which data are presented during the unsteady runs was 
generally within +0.2 foot per second (6.5 cm/s) and +0.2 revolution 


per second of the target values respectively. 
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For runs with fixed hull pitch angle w, (W=0), the value of p could 
be controlled to within +0.005 degree. For dynamic pitch runs +0, the 
selection of a propeller revolution at a specified ~ necessitated a toler- 
ance of 0.1 degree to ~; however, the average value of ~ for which data 
are presented during the unsteady runs was generally within 0.02 degree of 
the target }. 

Considering all sources of error including deviations during a run 
and inaccuracies in setting conditions, the model scale forces and moments 
presented in this report are generally considered to be accurate to within 


(plus or minus) the following variations: 


g Fax M MMAX 
1b = «© (N) 1b (N) jin-1b (N-m)}jin-1lb (N-m) 
Steady ahead V=0,U=0 | 0.1 (0.4) | 0.2 (0.9)| 0.2 (0.02)| 0.4 (0.06) 
Dynamic pitch V=0,040| 0.2 (0.9) | 0.3 (1.3)| 0.4 (0.04)| 0.6 (0.07) 
Acceleration V>0,=0 ONSH Gl3))E ON aaa Glee) Om One (ORO) a Ole (0.09) 


The values are somewhat more accurate for the steady-ahead runs than 
for the time-dependent runs, because the experimental conditions could be 
controlled more precisely for the steady runs and the measured forces and 
moments were averaged over many more revolutions of the propeller. The 
time-average values per revolution (based on 90 samples per revelolution) 
are slightly more accurate than the maximum values (based on one sample 
per revolution) which took into account the variation with blade angular 
position. Further, the peak values may have been slightly influenced by 
the dynamic response of the flexures, as discussed in the section on cali- 


bration. 


EXPERIMENTAL RESULTS 


LOADING COMPONENTS 

The basic loading components are shown in Figure 1. For a right- 
hand propeller the sign convention follows the conventional right-hand 
rule with a right-hand Cartesian coordinate system. For a left-hand pro- 
peller all the loads are the same, but for this case the sign convention 


follows a left-hand rule with a left-hand Cartesian coordinate system. 
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The sign conventions for both right-hand and left-hand propellers are 

shown in Figure 1. In all pertinent figures and tables throughout this 

report the blade loading components are listed in the following order: 
1. Components measured by Flexure 1: 

a. Ee - axial force, or thrust per blade. 

b. Me - bending moment about the axis normal to the shaft axis 
at r=0. This moment is generated primarily by the B component of 
force. 

2. Components measured by Flexure 2: 

a. By - tangential force, or force normal to the propeller axis 
and the spindle axis. 

b. M — moment about the propeller axis, (r=0), or torque per 
blade. This moment is generated primarily by the By component of 
force. 

3. Components measured by Flexure 3: 

a. Ee - radial force, or force parallel to the blade spindle 
axis. 

b. M, - moment about the spindle axis, or spindle torque. 

4. Supplemental components which were derived from the components 
listed above (not derived for all conditions). 


ae My 3 = Fe Gn -0.3R) cos 0.3 + By (x 
— bending moment applied on the spindle axis about the axis inter- 


Fy -0.3R) sin $5.3 


secting the spindle axis at r=0.3R and parallel to the expanded 


pitch line at r=0.3R. The My 3 vector as defined by the convention- 


al right-hand rule for a right-hand propeller (left-hand rule for 
left-hand propeller) intersects the plane normal to the propeller 


axis at the angle bo ane tan (P, 3/0-3 mD) and is directed so that 


a positive M 3 puts the face (pressure side) of the blade in 


0 


tension. 


be Ei SOSA EED -0.4R) cos do.4 4p a ee -0.4R) sin 5.4 


— bending moment applied on the spindle axis about the axis inter- 


ONE 


secting the spindle axis at r=0.4R and parallel to the expanded 


pitch line at r=0.4R. The M vector as defined by the conven- 


0.4 
tional right-hand rule for a right-hand propeller (left-hand rule 
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for the left-hand propeller) intersects the plane normal to the pro- 
peller axis at the angle $6 Tia tan (Py ,/ (0.47D)) and is directed 


so that a positive M puts the face (pressure side) of the blade 


0.4 


in tension. In calculating M from the experimental values of the 


0.4 
three measured forces and three measured moments, an adjustment was 
necessary to allow for the contribution of loading in the region 
0.4R>r>r, =0.3R where Th is the hub radius. It was estimated that for 
all harmonics including the time average values, 3 percent of M, and 
ah was contributed by the loading in the region 0.4R>r>r, These 
estimates were based on a refinement of the method of Cummings for 
the time-average values, and the method of Tsakonas et pile for the 
unsteady values. 

Hydrodynamic, centrifugal, and gravitational loads may contribute to 
each of these components of loading; however, for some components the 
centrifugal loads and/or gravitational loads are zero, as discussed in 
the section on centrifugal and gravitational loads. 

Each component of loading is generally presented as a variation of 


the instantaneous value with blade angular position ® and as a Fourier 


series in blade angular position in the following form: 


F,M(6) = (F,M) + ae (F,M) cos (né - ) 


(OF Mn 


where F,M = circumferential average value of F,M 
(F,M)_ = amplitude of the nth harmonic of F,M 
6 = angular position about the propeller axis, positive counter- 
clockwise from the vertical upward looking upstream for 
starboard propeller (left-hand rotation) positive clockwise 
looking upstream for port propeller (right-hand rotation) 
(OF wn = phase angle of nth harmonic of F,M 
where the reference line of the blade is the spindle axis; see Figure 2 
and Table 1. 


The components SS and M, are the most important for determination of 


the time-average and unsteady stresses in the hub mechanism of an actual 


BS) 


controllable pitch propeller. The components Mo 3 and My 4 are the most 
important for determination of the time-average and unsteady stresses in 


the blades of a propeller. 


CENTRIFUGAL AND GRAVITATIONAL LOADS 

The results of the air-spin experiments, corrected for interactions, 
are presented in Table 6. The time-average values arise from centrifugal 
force whereas the first harmonic arises from gravitational force. There- 
fore, the mean values should vary as ma where n is the propeller rotation- 
al speed, and the first harmonic should be independent of n. 

For the mean values, which arise from centrifugal force, signficant 
nonzero values were obtained only for the Hae Fo ee and M, components. 
Any realistic propeller would have nonzero values of centrifugal loading 
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components Go, and (Me: 


Nonzero values of oe and Coe are pro- 
duced by the nonzero values of skew and rake of the propeller evaluated. 
The components CO and mm), should be zero for any geometry, however, 
a small value of (Ey was measured. This small nonzero (Pade probably 
arises from inaccuracies in the air-spin experiment and interaction 
matrix. For all components the experimentally determined mean value 
varies essentially as i The experimental air-spin results were faired 
so that the values of the mean loading components used for separating 
hydrodynamic loads from total loads varied exactly as mee 

For the first harmonic loads, which arise from acceleration due to 
gravity, nonzero values were obtained only for the #09 M and Ee flexures. 
For Mie and Geile the phase angles are +96 degrees and -96 degrees, 
respectively; therefore the maximum and minimum values of these components 
occur when the blade is approximately horizontal. This would be expected 


from the geometry. The phase angle for ile is -159 degrees; therefore 


Ny seesTal R.J., "A Method of Calculating the Spindle Torque of a 
Controllable-Pitch Propeller at Design Conditions," David Taylor Model 
Basin Report 1529 (August 1961). 
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_ maximum value occurs when the blade is nearly vertical downward (6 
o'clock position). This is again as would be expected from geometry. The 
phase angles would not be expected to be precisely +90 degrees or 180 
degrees since the propeller has 22 degrees of skew. The amplitudes of 
Dig and Cag each should be equal to the combined weight of the 
blade and that portion of the appropriate flexure at radii greater than 
‘the radius of the appropriate strain gage. The values of Oia and 
ete wei confirmed by weighing the blade and appropriate flexures. 

The values of CO ag and Cae are essentially zero because the blade is 
skewed and raked so that mass of the blade is balanced about the spindle 
axis in both the plane containing the spindle axis and the propeller axis, 
and the plane normal to the propeller axis which contains the spindle 

axis (see Figure 2 and Table 1). The value of Cis is nearly zero since 
lies is always in a nearly horizontal direction. For all components the 
experimentally determined amplitude and phase of the first harmonic were 
essentially independent of rotational speed n. The experimental air-spin 
results were faired so that values of the amplitude and phase of the first 
harmonic of the loading components used for separating hydrodynamic loads 
from total loads were constant, independent of n. 

Approximate scaling parameters for centrifugal loads are (F/ppn°D") 
and (M/o,n°D>), whereas appropriate scaling parameters for gravitational 
loads are (F/opgD°) and (M/0p8D"). The model experiments presented in this 
report were conducted at full scale values of Froude number FA = (V/VveL) 
and advance coefficient J = (V/nD). Constant Froude number implies that 

V~(gL) ?~(gb)? 

V-~eD 
Constant advance coefficient implies that 

V~nD 

vee nee 
Therefore, 

ened 

ppeD’~ppn-D* 


5 
‘PpsD ~Ppn D 
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Thus, for the results presented in this report gravitational loads and 
centrifugal loads scale the same. Furthermore, if proper allowance is 
made for the difference in density between the model propeller and the 
full scale propeller,* the gravitational and centrifugal loads scale the 
same as the hydrodynamic loads. 

Therefore, in addition to the components of loading arising from 
hydrodynamic effects alone, for many experimental conditions the compo- 
nents of loading arising from the sum of hydrodynamic, centrifugal, and 
gravitational effects are presented. The components of loading arising 
from the sum of the hydrodynamic, centrifugal and gravitational effects 
are designated components of total loading. The centrifugal and gravi- 
tational loads presented are equivalent values for a nickel-aluminum- 
bronze propeller blade. These combined, or total, loading results are 
discussed in later sections. 


The time-average centrifugal spindle torque results, M, are com- 
c 
pared in Figure 9 with calculated values using the method of RELL 


Previous measurements of spindle torque by Boswell et anki and by Hawdon 
et ai have correlated well with values calculated by this procedure. 


The calculated value of M, is 33 percent lower than the experimental 
c 
value at design pitch (see Figure 9); however, this is within experimental 


accuracy. The largest measured value of M, is 0.5 inch-pounds (0.07 N-m) 
e 
as shown in Table 6 whereas the accuracy of this measurement is plus or 


minus 0.2 inch-pounds (0.02 N-m) as discussed in the section on accuracy. 


The large experimental inaccuracy as a percent of the measured M, value 
c 
results from a combination of (1) the small value of the measured M, ; 
c 
and (2) the inferior characteristics of flexure number 3, which measures 


*The model propeller used in these experiments was made of aluminum, den- 
sity pp=5.44 lbf-s2/£t4 (2.80 g/cm3). The full scale propeller on the 
DD-963 is made of nickel-aluminum-bronze, density Pp=14.48 lb£-s2/f£t4 
(7.46 g/cm3). 


re crieilil, R.J. et al, "Experimental Spindle Torque and Open-Water Per- 


formance of Two Skewed Controllable-Pitch Propellers," David Taylor Naval 
Ship Research and Development Center Report 4753 (December 1975). 
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EN and M,> relative to the other two flexures as discussed in the section 


on calibration. 


INFLUENCE OF DYNAMOMETER BOAT 

The results of the wake surveys with and without the downstream body 
(dynamomemter boat) are presented in Figure 10, and in Appendix A. These 
data indicate that the downstream body had only a small effect on the 
circumferential and radial variation in the flow and only a small effect 
on the harmonic content of the flow. However, they also indicate that the 
downstream body reduced the volume mean velocity through the propeller 
disk by approximately 12 percent; i.e., without the downstream body the 


volume mean wake (1-w)=1.06 and with the downstream body (1- m7 0-93- 


W. 
These results are, of course, without the propeller in place. i 

The change in effective velocity through the propeller due to the 
downstream body was deducted from thrust and torque identities between the 
mean thrust and torque measured during the blade loading experiments at 
the self propulsion point (Condition 1 in Table 3), and mean thrust and 
torque measured during a previous self propulsion model experiment.* 
These results, which include the effect of the propeller, indicate that 
the downstream body reduced the effective velocity through the propeller 
disk by approximately 5 percent; i.e., without the body, (1-w,)J=1.02 and 
Caton men sou: whereas, with the body, (1-w,,)=0.97 and Cheat Ooe2 

The difference between the effect of the downstream body on volume 
mean wake and effective wake is probably due to a combination of the 
following: 

1. The effect of the propeller action; (1-w,) and Cy) include 
the effect of the propeller but (1-wyy) does not. 

2. Experimental inaccuracies; both methods for calculating the 


change in velocity are based on a small difference of two much larger 


nearly equal experimental results. 


*DINSRDC experiments 21 and 22 on Model 5265-1B, in which the mean 
thrust and torque was measured using transmission dynamometers mounted 
inside the model hull. 
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It is judged that the dominant cause of the discrepancy is the effect of 
the propeller. 

Based on these results it was concluded that the downstream body 
reduced the mean velocity into the propeller by 5 percent at the self- 
propulsion condition. This is somewhat smaller than the 10 to 14 percent 
reduction in effective wake that was obtained in Reference 2 in which the 
same dynamometer boat was used behind a single screw model hull. It was 
assumed that the 5 percent reduction in the present experiment occurred 
at all conditions at which experiments were conducted. Therefore, after 
the effects of centrifugal force were subtracted from the measured loading 
components as discussed previously, the time-average value per revolution 
of each hydrodynamic loading component was corrected for the effect of the 
downstream body as follows: From the measured hydrodynamic blade thrust 


CE ) and hydrodynamic blade torque (M ), effective advance coefficients 
H 
based on thrust identity (Jp) and torque identity (J) were deduced from 


the open-water data (Figure 11). These values were multiplied by (1/0.95) 


to obtain corrected values of Jn and Jos i.e., without the downsteam body. 


The corrected values of La and M were then obtained from the open-water 
H 
data at the corrected advance coefficients Jr and Jae respectively. It 


was assumed that the downstream body did not affect the radial center of 


thrust F_ and tangential force F_. Therefore, 
Vy 
M corrected = (F_ corrected/F measured) (M_ measured) 
Vy aH Vy 


2 


-F corrected corrected/M measured) (F_ measured) 
Vy *y aH vy 


The spindle torque (M, ) was corrected by the same procedure as used 
for F_ and Gos using unpublished hydrodynamic spindle torque data for 


the DD-963 propeller. No corrections were made to Ee for the effect of 
H 
the downstream body; however, BD is very small for all experimental con- 
H 
ditions, as discussed later. 


No correction for the effect of the downstream dynamometer boat was 


made to the measured circumferential variation of the loading components. 
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Calculations made by the methods of Tsakonas et aie? and MeCArEnyaG indi- 
cated that the influence of the downstream body alters the peak-to-peak 
circumferential variation of the loads by no more than 2 percent. How- 
ever, these methods did not agree well with the experimental results, as 


discussed in the section on correlation with full-scale data and theory. 


STEADY-AHEAD OPERATION 

For operation near the self-propulsion point (Condition 1 in Table 
3), Figure 12 presents the variation of the various components of total 
blade loading with blade angular position and Figure 13 presents the 
amplitude of the first 25 harmonics of the various components of total 
blade loading. 

Based on the dynamic. calibration, as discussed in the section on 
calibration, it was judged that for all loading components the data are 
valid for the first 10 harmonics. In addition, the wake data shows no 
significant amplitudes for harmonics greater than the tenth; see Appendix 
A. Therefore, all data and analysis except Figures 12 and 13 are based 
on reconstructed signals using the first 10 harmonics. The symbols shown 
in Figure 12 indicate unfiltered values determined from the experiment; 
each represents the average value at the indicated blade angular position 
for over 200 propeller revolutions. The variation in measured values at 
a given angular position is discussed in the section on accuracy; see 
Figure 8. The lines on Figure 12 are the signals reconstructed from the 
first 10 harmonics. Figure 12 indicates that the variation of the 
signals with blade angular position are adequately represented by the 
number of harmonics retained. 

Figure 13 shows that there are no significant resonances for any of 
the loading components below the 23rd harmonic, which corresponds to 
(23)x(14.08)=324 hertz. This is higher than the lowest frequency reson- 
ance obtained in the results presented in Reference 2; i.e., 247 hertz. 
As discussed in the section on calibration, the higher frequency of the 
fundamental significant resonance obtained in the present experiment was 
anticipated because smaller and lighter blades were used in this experi- 


ment than were used in Reference 2. 
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The variation of all measured loading components with blade angular 
position for the self propulsion condition (Condition 1 in Table 3) is 
shown in Figures 14 and 15 for the hydrodynamic loads, and is shown in 
Figure 16 for the total (hydrodynamic, centrifugal and gravitational) 
loads. The amplitudes and phases of the harmonics of these loading compo- 
nents are presented in Figure 17 for the hydrodynamic loads and in Figure 
18 for the total loads. Appendix B presents tabulated values of all the 
data in Figures 14 through 18, and Table 7 presents a summary showing the 
maximum value, minimum value and first harmonic of each loading component. 
The values for each loading component are presented as decimal fractions 
of the time-average value of the corresponding loading component. These 
time-average values for both hydrodynamic loads and total loads are pre- 
sented in Table 8. 

These data show that for hydrodynamic loading the variation of all 
loading components was predominantly a once-per-revolution variation. The 


extreme values for all loading components, except F, and My» occurred near 


the angular position of the spingle, axis, 0=124 ek degrees; i.e., 34 
degrees beyond the horizontal. At these positions the blade tip is 
approximately 12 degrees beyond the horizontal. This suggests that the 
tangential component of the wake is the primary driving force; see Figure 
10. The extreme values of Fy and M, occur at up to 120 degrees after the 
extreme values of the other components. The reason for this variation in 
location of extreme values is not clear; however, it may be partially due 
to experimental inaccuracy with the FO-M, flexure as discussed in the 
section on calibration. 

For total loading, the variation of all components with blade 
angular position follows basically the same pattern as for hydrodynamic 
loading. This occurs because the unsteady loading due to gravity, which 
is a pure first harmonic of blade angular position, is much smaller than 
that due to hydrodynamic force. Further for all components with a mea- 
surable gravitational load, except Fo» the maximum value occurs near the 
angular position at which the spindle axis is horizontal; i.e., the 


gravitational load is nearly in phase or 180 degrees out of phase with 


the hydrodynamic load. 
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For hydrodynamic loading, F, and M_ were the largest measured force 
H H 
and moment components; see Table 8. For these components the maximum 


values were approximately 1.43 times the time-average values, and the 
maximum value minus the minimum value (double amplitude) was approximately 
0.91 times the time-average values; see Table 7. For F_ and|M |, the 

YH “A 
maximum values and range of values with blade angular position were 
slightly smaller fractions of the respective time-average values. For 
nae? the maximum value and range of values with blade angular position 
were much greater fractions of its time-average. This large fractional 


variation in F occurs because|F M Hee maximum 
Zz Zz Zz 
H H H 
value and range of values were 1.30 and 0.6/7, respectively, times the time- 


was very small. For 


average value. The radial point of application of FY varies from 0.68R 
H 
to 0.73R, and the radial point of application of Bo varies from 0.6/R to 
H 
O.79R. 


The maximum values of es and a were approximately 1.41 times the 
time-average values, and the range of values were approximately 0.88 times 
the time-average values for combined hydrodynamic, centrifugal and gravi- 
tational loading components, or total loading components. These results 
are nearly the same as the hydrodynamic results since the centrifugal and 
gravitational loads are small for these components; see Tables 6 and 8. 

The centrifugal and gravitational loads are a significant portion of 


the total loads for other loading components; for total loads F /EL=1.14 
MAX 
whereas for hydrodynamic loads F /F_=1.40. This large difference 
YuMAx YH 
results from the combination of centrifugal force adding to the time- 
average hydrodynamic force and gravitational force substracting from the 
circumferential variation of hydrodynamic force. Similarly, for total 


loads M /M,. =1.26 whereas for hydrodynamic loads M, /M 
X H,MAX “H 
centrifugal loading Ey is two orders of magnitude larger than the time- 


average hydrodynamic loading, and for M 


R634, ava 


7 the centrifugal loading is almost 


as large as the time-average hydrodynamic loading; see Table 8. 


The results presented here follow trends similar to those in Refer- 
ence 2 for the FF-1088 which is a single screw transom stern configura- 
tion. The component Oe which is the largest moment component for both 


cases, yields He ae =1.40 for the present configuration (DD-963 Class) 
MAX 
and M /M =1.38 for Reference 2 (FF-1088). The maximum and minimum 
YMAX yi) 
values occur at approximately the same angular position of the blade mid- 


chord at the 70 percent radius for the two configurations. 


HULL PITCH 

Figure 19 presents the variation of the peak values and time-average 
values per revolution of the various components of blade total (hydro- 
dynamic, centrifugal and gravitation) loading* with hull pitch angle yp for 
both quasi-steady simulation (time rate of change of hull pitch angle 
b=0) and unsteady simulation (W#0). These data show that for the quasi- 
steady simulation the time-average value per revolution of each loading 
component remains within 6 percent of its value corresponding to self- 
propulsion in calm water. The time-average value per revolution for the 
unsteady simulation, deviates by up to 12 percent from its value corre- 
sponding to self propulsion in calm water. 

Data at each specified value of hull pitch angle w for the quasi- 
steady runs were recorded and averaged for a minimum of 200 propeller 
revolutions whereas data for the dynamic pitching runs at each specified 
W represented an average of from 10 to 35 propeller revolutions. As dis- 
cussed earlier, the selection of a propeller revolution at a specified yp 
during the dynamic pitch runs necessitated a tolerance of only 0.05 degree 
to ~. Therefore, the differences between the results for the quasi-steady 
and unsteady simulations, including the time-average values per revolution, 
were Significantly larger than any errors which may have arisen from inac-— 
curacies in setting the experimental conditions. 

For quasi-steady simulation, the absolute value of the time-average 


value per revolution of all loading components, except spindle torque M,> 


+ 

No results are ahown for F_ since the F_ loading arises primarily from 
centrifugal effects, as discussed previousl¥, which are independent of 
hull pitch. 
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were larger for the stern-up condition than for the stern-down condition; 
the largest value occurs at (UP yp =1-85 degrees or time=0.31 seconds in 
the reference of Figure 19. This suggests that the effective speed of 
advance of the propeller increases slightly for the stern-down condition 
and decreases slightly for the stern-up condition. This appears reason- 
able since for stern-up the propeller tends to be further into the bound- 
ary layer of the hull. However, the time-average value per revolution did 
not monotonically increase with increasing ~ for all components. 

For dynamic simulation the largest absolute value of the time- 
average value per revolution of all loading components, except spindle 
torque My» occurs at approximately 0.15 second after the condition 
WV oy) =0> W>0, which is the reference for time t=0 in Figure 19. This 
indicates that the maximum time-average value during dynamic simulation 
occurs at a value of hull pitch angle ~ which occurs 0.16 second or 0.1 
cycle, before the » at which the maximum time-average value occurs during 
quasi-steady simulation. 

There was a significant difference between the peak values for the 
quasi-steady simulation and the unsteady simulation. For the quasi-steady 
simulation, the variation of the peak values with hull pitch angle yp 
followed approximately the same trends as the variation of time-average 
values per revolution. These quasi-steady results indicated that for 
V-Voy up to 1.85 degrees, the maximum increase in the peak value of any 
loading component above the corresponding value for vv oy was 5 percent. 
For the dynamic simulation, however, the maximum value of the peak loads 
increased as much as 23 percent above the corresponding value for steady 
ahead at a fixed hull pitch YVoye 

The dynamic simulation exhibited a dramatically different trend of 
peak load with ~ than was indicated by the quasi-steady simulation. For 
the dynamic simulation, the largest value of the peak loading, for all 
components except spindle torque My» occurred at approximately time t=0.8 
second in the hull pitch cycle shown in Figure 19. This corresponds to 
w=1.5 degrees stern down during the portion of the cycle in which the 


stern is moving down; i.e., (-Voy=-1-5 degrees, t<0. For dynamic 
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simulation, the smallest value of peak loading occurred near VV ony as the 
hull passed from the stern-down to the stern-up portion of the cycle; i.e., 
(W-W ay )=0, ¥>0. 

This difference in the unsteady loading between the quasi-steady and 
unsteady simulations may be due to an additional relative velocity compo- 
nent arising from the motion of the hull during dynamic pitching. As the 


hull passes through p=) the vertical velocity of the hull (and propel- 


> 
ler) is a maximum. As a hull goes from stern up to stern down through 
YVoyp the upward velocity component relative to the propeller in the 
plane of the propeller tends to increase above the values at fixed hull 
pitch at v=Voye This tends to increase the amplitude of the first harmo- 
nic of the tangential velocity, and thereby increase the unsteady loading 
(and increase the peak loading). The maximum vertical velocity of the 
propeller for sinusoidal pitching with Cbuax Yew 7d: 85 degrees and fre- 
quency=0.8 hertz is approximately 1.47 feet per second (0.448 m/s). This 


is equivalent to an additional tangential velocity ratio OND) Ose, (05 a13}3}., 


For ~ fixed at v=Voy ((V,),/V)=0.130 (see Appendix A). Therefore, 


(VM ax, i40 = 0.130 + 0.133 = 2.02% 
((V),/V) §-0,0-Voy 0.130 


This maximum occurs at a model simulated time of approximately 0.2 second 
before the maximum measured loads. The measured increase in unsteady loads 
arising from dynamic pitching was somewhat smaller than this calculated 


increase in tangential velocity, for example: 


F = 10 
. x ba 
MAX , 70 ca = LGD & i. /Mil 
: ay 0.4 


*A numerical error was found in a similar calculation presented in Ref- 
erence 2, With the numerical error corrected the results in Reference 2 
are substantially the same as those presented here. 
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On the basis of two-dimensional quasi-steady theory, the increase unsteady 
loading should be approximately proportional to the increase in tangential 
velocity.* 

The unsteady loading is important from consideration of fatigue of 
the propeller blades and hub mechanism. Since a ship may operate for an 
extended period in a seaway, the effect of the ship motions, such as 
dynamic hull pitching, on unsteady blade loads is significant. The differ- 
ence between the peak load and the time-average load per revolution is a 
measure of the unsteady loading. With this difference as a measure of the 
unsteady loading, the quasi-steady simulation indicates that for hull 
pitch angles b-Voy up to 1.85 degrees, the unsteady loading for Mo which 
is the largest moment component, increased by 8 percent above its corres- 
ponding value for v=Voye By contrast, the dynamic simulation showed the 
unsteady loading for the Mt component increased by 50 percent above its 
corresponding value for Vow without hull pitching. This indicates that 
the quasi-steady simulation is completely inadequate for estimating the 
effect of the seaway on unsteady loading. This also shows that the effect 
of the ship motions can dramatically increase the unsteady loading on the 
blades. Therefore, the effect of the ship motions due to operation in a 
seaway should be considered in any analysis of blade loading and in any 
fatigue analysis of the propeller blades or hub mechanism. 

The results presented here for hull pitching generally agree with 
the same type of results presented in Reference 2 for a model of the FF- 
1088, which is a single screw transom stern configuration. For ae which 
is the largest measured moment component in both cases, the comparative 
results, presented as a fraction of the time-average value without hull 


pitching, are as follows: 


“Chis simple analysis provides an upper bound to the dynamic pitching 
load, since the hull boundary above the propeller would tend to reduce 
the dynamic pitching induced upward tangential velocity relative to the 
propeller. 
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DD-963 FF-1088 


(Present (Reference 

Report) 2) 
Peak value, W#0 1.60 Lo S7 
Peak value, =0 aS 1.40 
Peak value without pitching v=0, v=Vow 1.40 1.36 
Maximum time-average value, wW#0 EO) 1.03 
Maximum time-average value, W=0 1.05 1.04 


The variation of the loading components with simulated time during the 
pitch cycle are also somewhat similar for these two configurations. The 
comparative results for Ho presented as time in seconds following the 


point YY oy=O> W>0 are as follows: 


DD-963 FF-1088 

(Present (Reference 

Report) 2) 
Peak value, W#0 0.77 0.62 
Peak value, W=0 0.31 0.31 
Maximum time-average values, #0 0.20 0.72 
Maximum time-average value, W=0 0.31 0.31 


ACCELERATION 

The variation of all measured loading components with blade angular 
position for the quasi-steady simulated acceleration condition Vv=P=n=0 
(Conditions 7 through 11 in Table 3) is shown in Figures 20 and 21 for the 
hydrodynamic loads, and is shown in Figure 22 for the total (hydrodynamic, 
centrifugal, and gravitational) loads. The amplitudes and phases of the 
harmonics of these loading components are presented in Figure 23 for 
the hydrodynamic loads and in Figure 24 for the total loads. Appendix B 
presents tabulated values of the data in Figures 20 through 24. The 
values for each loading component are presented as decimal fractions of 
the time-average value of the corresponding loading component at the self 
propulsion condition (Condition 1 in Table 3). These average values for 


both hydrodynamic loads and total loads are presented in Table 8. 
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Figure 25 presents the Taylor wake fraction based on thrust 1-w,, and 


ae 
the Taylor wake fraction based on torque rele as derived from the measured 
values of ie and M, and the open-water characteristics of the propeller 


H H 
(Figure 11). These data indicate that (1-w,,) varies by only approximately 


3 percent during the simulated acceleration. The value of Cy varies 
by only 1 percent for simulated time t>40 seconds; however, the value of 


(1l-w,) varies substantially during the initial portion of the simulated 


Q 
acceleration (t<40 seconds). 
Figures 20 and 22 show that for all measured hydrodynamic and total 


loading components, except Fy which is small, the peak values, including 


H 
variation with blade angular position occurred at the self propulsion con- 


dition. That is, for the acceleration condition simulated (see Figure 7 
and Table 3), the propeller is not exposed to higher peak loads than those 
to which it is exposed during full-power steady-ahead operation. 

In contrast to the peak loads, for most components the largest time- 
average loads per revolution occurred at the first experimental condition 
(V=2.65 knots, n=10.21 revolutions per second, J=0.64) during the simulated 
acceleration maneuver. The largest measured hydrodynamic force and moment 


components, F andM_, yield (F /E )=1.21 and (M /M = 
*y vy *y MAX “H,SP Yu,MAx YH,SP 

1.29, whereas for total loading (F pe )=1.16 and (M / Mou) =) leepoulie 

SP YMax Sp 

The conditions (FL /F ) > (F /F +) and (M /M ys 

ele pisie H,MAX “H,SP age YH,MAX H,SP 

(M /M ) occur because the centrifugal and hydrodynamic components are 

Ymax sp 

additive for Fy and Mt (i.e., they have the same signs) and the hydro- 

dynamic loads increase with decreasing rotational speed n, whereas the 

centrifugal loads decrease with decreasing n. 

Higher time-average and peak loads than those shown in Figures 20 
and 22 could, of course, be developed during acceleration maneuvers, 
depending on values of v, nh, and P. 

Figure 21 shows the variation in the radial center of longitudinal 


FOGCe is and radial center of tangential force, r . These results 


F 
iE Vy 


show that the time-average radial centers of these force components vary 


F 
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monotonically with advance coefficient over the range evaluated. As the 


advance coefficient based on thrust effective wake, Jp=Jy -wy) increases 


from 0.63 (at V=2.65 knots) to its design value of 1.14 (at V=6.52 knots), 


= decreases from 0.76R to 0.71R whereas i increases from 0.66R to 
mE Vy 
0.74R. This variation in he is the reason that (@ /E y< 
oh H,MAX “H,SP 
OL /M ) as discussed in the preceeding paragraph. 


Yu,MAX H,SP 

For all loading components, the variation with blade angular posi- 
tion tended to be dominated by the first harmonic for all conditions 
throughout the simulated acceleration maneuver. For all conditions at 
which there was significant variation in loading with blade angular posi- 
tion, the maximum and minimum values for all components except ES and Mm 
occurred for the blade spindle axis near 9=135 or 315 degrees (blade tip 
near §=115 or 195 degrees). This suggests that the variation in loading 
with blade angular position is produced primarily by the circumferential 
variation of the tangential velocity in the propeller plane (see Figure 
10). The angular variation of each loading component retained basically 
the same shape independent of speed and advance coefficient. 

There was a dramatic reduction in the circumferential variation of 
all measured loading components with decreasing speed V and decreasing 
rotational speed n. Previous data have shown that for a given propeller 
in a given flow field, the circumferential variation in the hydrodynamic 
loading varies approximately as the product of ship speed V and rotational 
speed n; see Wersidisma. = Figure 26 presents results in a form which 
allows evaluation of how closely the measured unsteady loading varies with 
nV. The ordinate is the first harmonic of the components of hydrodynamic 


blade loading except Fo ,» which is very small, and the abscissa is nV. 
H 
The data shown in Figure 26 indicate that the first harmonic of each of 


the presented hydrodynamic loading components is approximately proportion- 


al to nV. 


= isseeiidieme R., "Tendencies of Marine Propeller Shaft Excitation," 
International Shipbuilding Progress, Vol. 19, No. 218, pp 328-332 (Octo- 
ber 1972). 
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Figure 27 presents the variation of the time-average values per 
revolution and peak values of the various components of total blade 


loading for both quasi-steady simulated acceleration (V=n=P=0) and unsteady 
simulated acceleration (V>0, n=P=0). 

There was only a small variation in the measured loading components 
between the quasi-steady simulated acceleration and the unsteady simulated 
acceleration. For all components except M,» the largest variation between 
the results from the two types of simulation expressed as a decimal fraction 
of the corresponding time-average value at the self-propulsion point was 
0.05 for the peak values and 0.02 for the time-average value per revolu- 


tion. The corresponding variations for M, were no greater than 0.06 for 


the peak values and 0.05 for the time-average value per revolution. 

The variation in the results between the two types of simulation 
appeared to be essentially random. This suggests that these deviations 
are some measure of the experimental accuracy and do not represent any 
systematic trends arising from the difference in V between the two types 
of simulation. 

Data for the quasi-steady simulation were recorded and averaged for 
a minimum of 200 propeller revolutions, whereas data presented for the 
unsteady runs represent an average of only five revolutions. Further, 
the steady experimental conditions which were set during the quasi-steady 
simulation allow the values of V and n to be controlled more precisely 
than during the unsteady runs; however, the average of the five values of 
V and n during the unsteady runs for which data are presented was general- 
ly within one percent of the target values. 

The results presented in this section for a simulated acceleration 
maneuver follow trends similar to those in Reference 2 for a simulated 
crash~forward maneuver on a model of the FF-1088. Both sets of data show 
the following: 


1. The values of (1-w,) and (l-w.) do not vary substantially except 


Q 
during the initial stages of the acceleration or crash-forward maneuver. 
2. The variation of all loading components with blade angular posi- 


tion was dominated by the first harmonic throughout the simulated maneuver. 
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3. The amplitude of the first harmonic of all loading components 
varied essentially as nV. 

4. The acceleration of the hull did not have a significant effect 
on the loads; i.e., the loads for quasi-steady acceleration V=0 and un- 
steady acceleration v>0 were not significantly different. 

The ratics of either the peak or time-average loads during the accelera- 
tion (or crash-forward) maneuver to the time-average loads at the self- 
propulsion point do not agree closely for the results in the present 
report (DD-963 Class) and Reference 2 (FF-1088). This difference is to 
be expected since these ratios are very sensitive to the value of V, n, 
and P for the simulated maneuvers, which are quite different for these 


two cases. The largest moment component for both cases ae gives 


/M = 1.40 from the present report and 1.51 from Reference 2. The 
Ypeak Ygp 
results from the present report show M /M = 1.21 whereas the results 
y. My, 
a i MAx ‘sp 
from Reference 2 yield M /M SValsos 
Ymax Ysp 


CORRELATION WITH FULL-SCALE DATA AND THEORY 

For operation near the self-propulsion point (Condition 1 in Table 
3), correlation was made between the model experimental loads obtained in 
the present investigation, bending moments deduced from strains measured 
on the corresponding full-scale propeller, and analytical calculations. 

These comparisons were made for My 3 and Moo 4 which were calculated 
from the three measured forces and three measured moments. As discussed 
in the section on experimental results, My .3 is defined as the bending 
Moment applied on the spindle axis about the axis intersecting the spindle 
axis at r=0.3R and parallel to the expanded pitch line at r=0.3R. The 
My 3 vector as defined by the conventional right-hand rule for a right- 
hand propeller (left-hand rule for a left-hand propeller) intersects the 
plane normal to the propeller axis at the angle 9 3 = tan "(Py 4/(0.3nD)) 


and is directed so that a positive M 3 puts the face (pressure side) of 


0 


the blade in tension. The component M is defined in a manner analogous 


0.4 


to My 3 except it, is referred to r=0.4R. 
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The full-scale data used for correlation are preliminary values of 
strains ¢« measured* at several chordwise stations at r=0.35R and r=0.45R 
on both sides of the blade on the DD-963 CP propeller. Both time-average 
strains per revolution and variation of strain with blade angular posi- 
tion were recorded during full-power stready-ahead operation. By inter- 
polation, values of radial strain at the spindle axis at r=0.4R were ob- 
tained. Assuming that the variation in radial strain is proportional to 
the variation in total (hydrodynamic, centrifugal, and gravitational) 


bending moment; i.e., that (e /[€ ) = /Mo 4)? these full- 


(M 
10.4 ax r0.4 0. 4vrax 


scale data indicate that (M /M. ,) = 1.48. From the model data 
0. 4uax 0.4 


[Me -)— ll W43and) OM 
0.4 say 


( 0.3 


M /M,. ,) = 1.56. A cursory evaluation 
0. 3ax 0.4 


of the variation of the full scale strain with blade angular position indi- 

cates that it follows trends similar to the bending moment determined from 

the model experiments. The correlation with full scale data presented 

here is preliminary; a more thorough correlation with the full scale data 

will be undertaken when analysis of the full scale data is complete. 
Theoretical calculations of hydrodynamic loads were made by using 


the method of Tsakonas et afi, 22 


which is based on unsteady lifting- 
surface theory and the method of Mecarehy co a quasi-steady technique 
which utilizes the open-water characteristics of the propeller. Although 
the method of McCarthy is a simple quasi-steady technique, it was judged 
that this method should be suitable to the cases under consideration in 
this report because the dominant unsteady loading occurs at a low reduced 


frequency and the dominant first harmonic of the wake is in phase radially 


*The full-scale measurements were conducted by personnel in DTINSRDC 
Code 1962 under the direction of G.P. Antonides. The results presented 
here are preliminary, and thorough analysis of the data is continuing. 
The details of this full-scale trial will be reported in a future DINSRDC 
report. 


a ieaoncice S. et al, "Documentation of a Computer Program for the 
Pressure Distribution, Forces and Moments on Ship Propellers in Hull 
Wakes," (In Four Volumes), Stevens Institute of Technology, Davidson 
Laboratory Report SIT-DL-76-1863 (January 1976). Revised April 1977. 
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(see Figure 10 and Appendix A). These calculations were made for Condi- 
tion 1 in Table 3 using the wake measured in the plane of the propeller 
both with and without the downstream dynamometer boat in place (Figure 10 
and Appendix A), and with the mean velocity through the propeller deter- 
mined from thrust identity used as the reference velocity. 

The use of the speed of advance based on thrust effective wake, 
V,=V(1-w,) as the reference speed in these calculations is consistent 
with the use of this velocity to correct the time-average loads for the 
effect of the dynamometer boat as discussed in the section on experimental 
results. Tsakonas et Bae recommend using the ship speed as the refer- 
ence velocity, which is equivalent to using (1-w,,)=1.05 however, this 
recommendation was not followed here because the flow does not pass through 
the propeller at the ship speed. To evaluate the sensitivity of the pro- 


2 : 
252) to the reference speed, calculations were 


cedure of Tsaknoas et al, 
performed for the first harmonic using the thrust effective wake (1-w,p) 
and the volume mean wake determined from the wake surveys (1-wy)) + These 
calculations showed the following: 


Wake without dynamometer boat 


M (using (1-w,,) = 1.02) 
OPS rene oie sal = 0.99 

HO-3e (using (1-w,,,) = 1.06) 

O04 by *yo.3, Ou = -0.3 degrees 

Wake with dynamometer boat 
M (using (1-w,,) = 0.97) 
si ie = 1.02 
HO. 34 (using (1-w,) = 0.93) 


SMOuse nT a *yo.3, uw? = 0.5 degrees 


Therefore, the calculated unsteady loads using the method of Tsakonas et 
25,34 
Allain are not sensitive to the reference speed over the range of con- 


cern in the present case. 
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22S calculations were conducted 


For the method of Tsakonas et al, 
for the first ten harmonics of the Tae The "normal'' components of wake 
harmonics, as required by this method, were defined as the wake harmonics 
normal to the chord line of the blade section at the local radius rather 
than normal to the advance angle at the local radius as recommended by 


oh) With the wake harmonics resolved normal to the blade 


Tsakonas et al. 
chord, this method apparently considers both the unsteady flow normal to 
the resultant inflow and an approximation to the unsteady flow parallel to 
the resultant inflow. 

The quasi-steady calculations are based on the circumferential vari- 
ation of the wakes measured at the 0.7 radial station. These calculations 
were made at 10-degree increments of blade angular position ®. It is 
assumed that the radial centers of the unsteady thrust and tangential 
force are at r/R=0.70 for all blade angular positions. 

Figure 28 presents values of HO. 3y and one with and without the 
downstream dynamometer boat, calculated with the methods of Tsakonas et 
aio and Me Carew, Based on these calculated results it appears that 
the dynamometer boat does not have a significant influence on the cir- 
cumferential variation of the blade loads. 


Figures 29 and 30 present the variation with blade angular position 


and the first ten harmonics of My 3 and My I from the model experiments 
34 “4 


and analytical calculations. All data are nondimensionalized on the same 
quantity, i.e., the time-average bending moment determined from the model 
experiments and corrected for the downstream body as discussed in the sec-— 
tion on experimental results. This comparison indicated that the experi- 
mental unsteady hydrodynamic bending moments were substantially higher than 


the calculated results. A typical comparison is as follows: 


*These calculations were made by using the computer program developed by 
the Davidson Laboratory including refinements made through April 1977. 
None of the refinements made since December 1975 influence the calculated 
unsteady loads presented in this report. Therefore, this calculation pro— 
cedure is the same as that used in Reference 2 for calculating the unsteady 
bending moments on the propeller on the FF-1088. 


Tsakonas, S. et al, "Correlation and Application of an Unsteady Flow 
Theory for Propeller Forces," Transactions of the Society of Naval Archi- 
tects and Marine Engineers, Vol. 75, pp 158-193 (1967). 
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Prediction Method 


M -M M —-M 
0 34 MAX 0.3, 0.45 MAX 0 4 
My 3 (experiment) My 4 (experiment) 

H H 
Model Experiment 0.36 0.36 
Quasi-Steady 
Procedure26 0.25 0.24 
Unsteady Pregecdeae = 0.20 0.19 


For this typical comparison the experimental result is approximately 47 
percent higher than the quasi-steady prediction and approximately 85 per- 
cent higher than the unsteady prediction. 

The circumferential variations in the model experimental results of 


other components of blade loading F_ , F_ ,M _,andM_ were larger than 
x y x ny; 
H H H H 

the values calculated by the two indicated procedures by approximately 


the same ratio as shown for My 3 and Mo 4 . These comparisons are not 
H H 
shown. 

Previous investigators have compared experimental unsteady forces 
and moments on a single blade of various propellers in inclined flow with 
forces and moments calculated by a quasi-steady procedure similar to that 
described by MEGAE Uy Mae These, experimental loads were obtained by direct 
measurement of unsteady forces and moments on a single blade (References 
2, 3, 17, 18, and 19) or were deduced from measured steady transverse 
forces and moments along axes fixed relative to the flow, i.e., not rotat— 
ing with the propeller (Reference 36). References 2, 3, 17, 18, 19, and 
36 all show that for noncavitating conditions, the experimental unsteady 
blade loading was from 1.5 to 2.0 times as large as the values calculated 


by the quasi-steady method. This agrees with the results of the present 


investigation; see Figures 29 and 30. 


eo enrachen F., "The Study of Ships' Propellers in Oblique Flow," 
Defence Research Information Centre Translation No. 4306, Copyright Con- 
troller: Her Majesties Stationary Office, London, England, October 1975; 
English Translation of "Untersuchung von Schiffsschrauben in schrager 
Anstromung,"’ Schiffbauforschung, Vol. 3, No. 3/4, pp 97-102 (1964). 
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Preliminary results from blade loading oaecimenne " conducted in 
idealized flows indicate that: 

1. The unsteady blade loads in either axial or tangential wakes are 
not influenced by the presence of a nearby boundary above the propeller. 

2. In inclined flow without an upstream hull, the experimental un- 
steady loading near the design advance coefficient is nearly two times as 


25,34 


large as that calculated by the method of Tsakonas et al, and approxi- 


mately 80 percent larger than that calculated by the method of MeCaceny oo 
3. In an axial wake with a dominant first harmonic of blade angular 
position which was generated by upstream wire grid screens with zero shaft 
angle, the unsteady loading near the design advance coefficient is within 
approximately 15 percent of the values calculated by the methods of 


25,34 


2 
Tsakonas et al, and McCarthy. ° This is in agreement with previous 


correlations with the method of Tsakonas et al for unsteady shaft (bearing) 
forces and moments for operation in axial calass Le 
These results indicate that the large discrepancy between the experi- 
mental and calculated unsteady bending moments presented in the current 
report appear to be due to the inability of the present theories to 
properly account for all important characteristics of the flow field for 
operation in inclined flow. All available procedures, including the un- 


2252) and the quasi-steady procedure of 


steady theory of Tsakonas et al 
vieGaawy implicitly assumed that the propeller slipstream follows the 
propeller axis rather than the direction of the effective velocity into 
the propeller which is at an angle to the propeller axis in inclined f.iow. 
It is speculated that this failure to consider the true direction of the 


slipstream is the major factor in the analytical underprediction of the 


*The results in Reference 12 were at substantially higher reduced fre- 
quency than the results in the present study; therefore, the quasi-steady 
procedure of McCarthy over-predicted the unsteady loads in Reference 12. 


Boswell, R.J. and S.D. Jessup, "Experimental Determination of Period- 
ic Propeller Blade Loads in a Towing Tank,'' Presented to the 18th American 
Towing Tank Conference, U.S. Naval Academy, Annapolis, Maryland (August 
LOTT) o 
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unsteady blade loads in inclined flow. Numerical computations to check 


this hypothesis are planned. 


SUMMARY AND CONCLUSIONS 

Experiments were described in which the mean and unsteady loads, 
including hydrodynamic, centrifugal, and gravitational loads, were mea- 
sured on a single blade of a model of a CP propeller on the DD-963 Class 
Destroyer. The experiments were conducted behind a model of the DD-963 
hull under steady-ahead operation, hull pitching motions, and simulated 
acceleration maneuvers. The discussion of experimental techniques in- 
cludes a description of the dynamometer and data analysis system. The 
results are summarized as follows: 

1. For all significant loading components, except for radial force, 
the loads are predominantly of hydrodynamic origin. 

2. The circumferential variations of all measured components of 
hydrodynamic and total blade loading are primarily a first harmonic, with 
maximum and minimum values occurring near the blade angular position which 
is 25 degrees past the position at which a radial line from the propeller 
axis to the tip is horizontal. 

3. For steady-ahead operation: 

a. The maximum values and peak-to-peak circumferential varia- 
tions for measured hydrodynamic forces and bending moments were up 
to approximately 1.43 and 0.91 of the time-average values, 
respectively. 

b. The maximum values and the peak-to-peak circumferential 
variations for measured total forces and bending moments were up to 
approximately 1.41 and 0.88 of the time-average values, respectively. 

c. The model results for circumferential variation of bending 
moments about the nose-tail lines of the 0.3 and 0.4 radii agreed 
fairly well with loads deduced from strain measurements on the full- 
scale propeller, but they were larger than theoretically calculated 


values. 
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4. For simulated hull pitching (maximum pitch angle of 1.85 
degrees): 

a. The maximum values of measured total forces and bending 
moments increased over the corresponding values without hull pitch 
by 5 percent for quasi-steady simulation and by 23 percent for un- 
steady simulation with model pitching frequency equal to 0.8 hertz 
(full scale equivalent frequency is 0.16 hertz). 

b. The peak-to-peak circumferential variation of the measured 
total forces and bending moments increased over the corresponding 
values without hull pitch by approximately 5 percent for quasi- 
steady simulation and by approximately 50 percent for unsteady 
simulation with model pitching frequency equal to 0.8 hertz. MThere- 
fore, any quasi-steady simulation of ship motions is completely in- 
adequate for estimating the effect of ship motions on unsteady pro- 
peller blade loading. 

5. For the simulated acceleration maneuver: 

a. The dominant first harmonic of the measured hydrodynamic 
forces and bending moments varied in a nearly linear manner with 
the product of ship speed and propeller rotational speed. 

b. The acceleration of the hull did not have a significant ef- 
fect on the measured loads. Therefore, propeller blade loading 
during an acceleration maneuver can be adequately estimated by 
quasi-steady experiments. 

c. The maximum time-average values of measured forces and bend- 
ing moments per revolution were in the range of 1.21 to 1.29 of the 
time-average values during full-power steady-ahead operation for 
hydrodynamic loads, and in the range 1.16 to 1.21 of the time-average 
values during full-power steady-ahead operation for total loads. 

d. The simulated acceleration condition did not expose the pro- 
peller to higher peak loads than those to which it is exposed during 
full power steady-ahead operation. However, these loads are very 
sensitive to the maneuver simulated and substantially higher peak 


loads could be developed during other acceleration maneuvers. 
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e. Except for the initial portion of the simulated accelera- 
tion maneuver, the Taylor wake fractions were within three percent 
of their values at the self propulsion point. 

All of the results presented here on a model of the DD-963 Class 
Destroyer follow close to previously reported results of similar experi- 


ments on a model of the FF-1088. 
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APPENDIX A 


DETAILS OF WAKES* 


Tables 9 and 10 present the velocity component ratios and the 
harmonic content of the wakes in the plane of the propeller, both with 
and without the downstream dynamometer boat. The data at even radial sta- 
tions were obtained by interpolation and extrapolation of the measured 


data as described in-Reference 38. 


*Al11 data presented in Appendix A were obtained from wake surveys con- 
ducted by R.F. Roddy, DINSRDC Code 1524. Further details of these wake 
surveys will be presented in a future DINSRDC report. 

2S oiyeae H.M., "Analysis of Wake Survey of Ship Models - Computer 


Program AML Problem No. 840-219F,"' David Taylor Model Basin Report 1804, 
March 1964. 
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APPENDIX B 


DETAILED EXPERIMENTAL RESULTS 


Table 11 presents detailed experimental results, including variation 
with blade angular position and harmonic analyses, for steady-ahead opera- 
tion at V=6.52 knots, n=14.08 rev/sec. The data in Table 11 are tabulated 
values of the data presented in Figures 12, 13, 14, 16, 17, and 18. 

Tables 12 to 15 present detailed experimental results including 
variation with blade angular position and harmonic analyses, for the quasi- 
steady acceleration conditions. The data in Tables 12 to 15 are tabulated 


values of the data presented in Figures 20, 22, 23, and 24. 
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Figure 5 - Typical Strain-Gaged Flexure 
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Figure 10 (Continued) 
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Figure 12 - Influence of Extraneous Signals on Measured Loads 
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Figure 12 (Continued) 
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Figure 12 (Continued) 
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Figure 12 (Continued) 
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Figure 12 (Continued) 
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Figure 13 - Experimental Data Showing Extraneous Higher Harmonics 
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Figure 17 - Harmonic Content of Experimental Hydrodynamic 


Loads for Steady—-Ahead Operation 
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Figure 20 - Variation of Experimental Hydrodynamic Loads with 
Angular Position for Quasi-Steady Acceleration 
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Figure 20 (Continued) 
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Figure 20 (Continued) 
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Figure 20 (Continued) 
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Figure 20 (Continued) 
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Figure 22 - Variation of Experimental Total Loads with Angular 
Position for Quasi-Steady Acceleration 


5 KNOTS 

.21 REV/SEC 
3.55 KNOTS 
10.96 REV/SEC 


6 
2 


ee ae 


IS 
V = 5.36 KNOTS \ a 
n = 12.70 REV/SEC S 
V = 6.26 KNOTS 
n = 13.78 REV/SEC eee 


V= 6.52 KNOTS 
n = 14.08 REV/SEC 


135 180 225 270 315 360 


0 45 90 
POSITION ANGLE 6 (DEGREES) 


Figure 22a - F 
x 


109 


My/Myop 


1.8 


1.6 


1.4 


0.8 


0.6 


0.4 


0.2 


Figure 22 (Continued) 


oO 


V= 2.65 KNOTS 
n = 10.21 REV/SEC 


V= 3.55 KNOTS 
n = 10.96 REV/SEC 


n = 13.78 REV/SEC 


| 


V= 6.52 KNOTS 
n = 14.08 REV/SEC 


i) 
: J 
\SS By 
V= 5.36 KNOTS 
n = 12.70 REV/SEC : | 
V= 6.26 KNOTS a 
4 


90 


135 180 225 
POSITION ANGLE @ (DEGREES) 


Figure 22b - Hh 


110 


270 315 360 


Figure 22 (Continued) 


= 6.52 KNOTS 
n = 14.08 REV/SEC 


ine 


V= 6.26 KNOTS 
n = 13.78 REV/SEC 


mo 


5.36 KNOTS 
n = 12.70 REV/SEC 


Fy/Fyop 


V= 3.55 KNOTS 
n = 10.96 REV/SEC 


10.21 REV/SEC 


2.65 KNOTS | 


0 45 90 135 180 225 270 315 


POSITION ANGLE 6 (DEGREES) 


Figure 22c - Bo 


111 


360 
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Figure 22 (Continued) 
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Figure 22 (Continued) 
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Figure 23 - Harmonic Content of Experimental Hydrodynamic 
Loads for Quasi-Steady Acceleration 
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Figure 26 - Variation of First Harmonic of Experimental Hydrodynamic 
Loads with nV for Quasi-Steady Acceleration 
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Figure 27 (Continued) 
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Figure 27 (Continued) 
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Figure 28 - Variation of Hydrodynamic Bending Moment at 30 Percent 
and 40 Percent Radii with Blade Angular Position, Theoretical 
Prediction With and Without Dynamometer Boat 
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Figure 28 (Continued) 
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Figure 29 - Variation of Hydrodynamic Bending Moment at 30 Percent 
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Comparison of Model Data with Theory 
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Figure 30 - Harmonic Content of Hydrodynamic Bending Moment at 30 Percent 
and 40 Percent Radii, Comparison of Model Data and Theory 
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Figure 30 (Continued) 
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TABLE 1 - CHARACTERISTICS OF PROPELLERS ON DD-963 CLASS 


Diameter, D: 17.0 feet (5.577 m)* 


Rotation: Inward Turning* 


Number of Blades, Z: 5 


Maximum Rotational Speed (Rated), n: 


r/R 


0.3 

0.35 
0.45 
0.55 
0.65 
0.75 
0.85 
0.90 
0.95 
1.0 


* For model propeller, D = 0.6848 feet (22.47 cm) 


168 r/min 


Full Power (Rated), Pp: 
40,000 hp (29,800 kW) 


c/D 


0.178 
0.210 
0.271 
0.327 
0.374 
0.406 
0.409 
0.387 
0.326 
0 


P/D 


1.165 
1.296 
1.480 
1.566 
1.566 
1.498 
1.381 
1.306 
1.222 
1.128 


0 ** (deg) 


2.985 
3.481 
4.810 
6.631 
8.978 
11.895 
15.410 
17.403 
19.557 
21.876 


DESTROYER; DINSRDC MODEL PROPELLERS 4660 AND 4661 


Speed at Full Power, V: 32.5 knots 
Hub-Diameter Ratio, D, /D: 0.30 
Expanded Area Ratio: 0.73 

Blade Thickness Fraction: 0.054 
Section Meanline: NACA a=0.8 


(Z,/D)** t/D fy/e 
-0.0006 0.0420 0 
-0.0022 0.0372 0.0050 
-0.0095 0.0290 0.0209 
-0.0185 0.0226 0.0267 
-0.0288 0.0178 0.0256 
-0.0392 0.0146 0.0209 
-0.0488 0.0122 0.0151 
-0.0529 0.0110 0.0122 
-0.0561 0.0091 0.0094 
-0.0583 0 0 


*The blade loads were measured on the starboard propeller (left-hand rotation); 


DTNSRDC Model Propeller 4661 


**The spindle axis is the blade reference line. 
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TABLE 2 — CALIBRATION MATRIX 


-3.33 -0.099 0.058 -0.002 -0.090 -0.066 
-0.008 5.00 -0.002 0.084 -0.168 -0.196 
-0.071 0.048 -3.34 0.070 -0.014 0.058 
Calibration Matrix = (C; ,) = 
J 0.009 -0.033 0.103 4.99 -0.071 0.297 
-0.234 -0.367 0.291 -0.861 -9.99 -1.87 


0.079 -0.103 0.223 0.020 -0.034 -5.02 


where, 
Fy Fx Fy Fy F, are indicated forces in volts 
My, My, My My. M,, are indicated moments in volts 
y) 2 Fyn sake Boe Fyn FA are applied forces in pounds 
Mx, Mx a vd Mixa! Mya! Me, are applied moments in inch- pound 
2) Fan C; for j= 1, 3,5 are in volts per pound 
M_, Ma C;; for j = 2, 4,6 are in volts per inch-pound 


The arrays are arranged by the coupling of the force and moment pairs of the three flexures; 
i.e., Flexure 1 measures F, and My, Flexure 2 measures Fy and M,, and Flexure 3 measures 
F, and M,. 
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TABLE 3 - MODEL EXPERIMENTAL CONDITIONS 


Condition Mie n V-Vow Vv 
Number | knot r/s (P/D)q 5 deste | knot/s 


Quasi-Steady 6.52 (3.61) 
Hull Pitch 6.52 (3.61) 
6.52 (3.61) 
6.52 (3.61) 


Quasi-Steady 
Acceleration 


Unsteady 2.65** (1.47) ** 
Acceleration 3.55** (1.96) ** 


(0.31) ** 
(0.31)** 
(0.18)** 
(0.04) ** 
(OES 


536 ne 2:97/) na 
6.26** (3.47) ** 
6.52** (3.61) ** 


oo o ooiooo0ocecodco 


*Sinusoidal with amplitude equal to 1.85 degrees, frequency equal to 0.8 Hz. 


**Varies with time (Figure 7); value shown is at time of interest. 


TABLE 4 —- FULL-SCALE CONDITIONS SIMULATED BY MODEL EXPERIMENTS 


Condition Wow 
cee knot apt oY ve patio 


Sat Proper 
Hull Pitch 


Quasi-Steady 
Acceleration 


(0.31) ** 
(0.31) ** 
(0.18) ** 
(0.04) ** 
(0) ** 


Unsteady 
Acceleration 


0 
0 
0 
0 
0 
0 
0 
0 
0 
0 


*Sinusoidal with amplitude equal to 1.85 degrees, frequency equal to 0.16 Hz. 


**\/aries with time (Figure 7); value shown is at time of interest. 
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TABLE 5 — REPEAT RUNS FOR BY FOR STEADY-AHEAD OPERATION 


Date of Run: 
Time of Run: 1100 
Run No.: 45 


Propeller Speed, 14.08 
n (r/s) 


Velocity, V 11.04 
(ft/s) 
n 
(1) 4.380 
1 1.628 
2 0.306 
3 0.215 
4 0.115 
5 0.060 
6 0.028 
7 0.008 
8 0.015 
9 0.009 
10 0.008 
1 0.008 
12 0.013 
13 0.009 
14 0.006 
15 0.005 
16 0.007 
1 118.8 
2 12.0 
3 38.2 
4 37.7 
5 31.9 
6 49.7 
7 167.9 
8 -126.4 
9 -135.2 
10 172.9 
11 155.3 
12 141.7 
13 133.7 
14 116.9 
15 90.6 
16 77.9 


18 Dec 76 


1300 


118.6 
11.2 
37.3 
36.8 
30.4 
46.5 

- 169.0 
- 125.5 
-126.3 
-178.3 

161.8 

140.9 

135.6 

112.7 
98.8 
78.4 


_ “Run 115 used for detailed analysis 
**Raw data without corrections for interactions or downstream body 


1410 
55 


14.08 


11.03 


4.408 
1.632 
0.305 
0.211 
0.117 
0.059 
0.024 
0.008 
0.017 
0.011 
0.009 
0.009 
0.014 
0.010 
0.007 
0.006 
0.008 


119.0 
11.6 
37.9 
36.6 
30.4 
39.8 

- 166.5 
-121.0 
- 133.4 
- 176.1 

176.2 

143.8 

143.2 

117.6 

113.7 
88.6 


1740 
57 


14.08 


11.07 


4.520 
1.532 
0.340 
0.221 
0.110 
0.054 
0.028 
0.007 
0.008 
0.010 
0.009 
0.012 
0.015 
0.012 
0.008 
0.003 
0.005 


115.3 
3.5 
19.1 
17.6 
-0.6 
46 
107.7 
-171.6 
140.1 
107.8 
106.7 
91.5 
95.4 
75.5 
110.8 
55.1 


21 Dec 76 
1052 1613 1120 
65 92 114 
14.08 14.01 14.08 
11.05 11.04 10.97 
Amplitudes (F)** (Ib) 
4.635 4.528 4.321 
1.631 1.622 1.563 
0.335 0.330 0.255 
0.228 0.234 0.201 
0.118 0.123 0.119 
0.063 0.069 0.054 
0.031 0.033 0.021 
0.003 0.002 0.003 
0.010 0.008 0.010 
0.002 0.004 0.003 
0.010 0.011 0.010 
0.011 0.011 0.012 
0.014 0.015 0.017 
0.010 0.011 0.012 
0.007 0.007 0.008 
0.002 0.004 0.005 
0.003 0.003 0.005 
Phases (b,)** (deg) 
116.8 116.5 116.7 
-0.0 -0.2 1.2 
24.1 23.8 21.6 
20.9 19.2 24.8 
10.5 9.5 13.5 
13.5 16.4 11.2 
113.7 116.5 | -159.1 
-155.3 179.3 | -137.8 
171.2 136.1 -157.9 
128.3 125.0 130.7 
130.6 111.9 127.5 
112.1 100.1 114.5 
114.5 108.3 123.8 
113.2 97.0 84.3 
96.1 74.4 60.3 
74.6 61.0 83.8 
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28 Dec 76 —————___ | 30 Dec 76 
1310 1315 1335 1000 
115° 116 117 128 


14.08 14.02 14.08 14.00 


10.97 10.97 10.97 10.97 


4.495 4434 4457 | 4.528 
1.608 1.594 1.588 | 1.610 
0.257. 0.261 0.255 | 0.292 
0.203 0.204 0.201 0.229 
0.123 0.123 0.121 0.127 
0.055 0.051 0.052 | 0.056 
0.025 0.022 0.021 0.027 
0.002 0.002 0.003 | 0.005 
0.010 0.011 0.011 0.013 
0.003 0.004 0.005 | 0.008 
0.009 0.009 0.010 | 0.010 
0.014 0.014 0.016 | 0.012 
0.016 0.017 0.019 | 0.018 
0.011 0.011 0.014 | 0.014 
0.006 0.005 0.008 | 0.008 
0.006 0.005 0.005 | 0.006 
0.005 0.006 0.006 |! 0.006 


116.5 116.3 116.2 118.7 
1.3 3.4 3.2 -2.7 
21.4 21.7 20.7 26.0 
23.6 26.5 23.7 30.3 
12.6 15.5 12.2 18.7 
5.8 6.2 3.0 35.5 
-92.1 161.0 -168.8 101.7 
-123.6 -127.9 -140.6 | -147.2 
-113.9 -120.8 127.9 172.0 
123.6 119.1 120.3 143.5 
118.4 122.3 121.3 146.1 
115.0 114.3 122.1 143.7 
121.2 119.7 123.2 132.1 
78.3 85.8 85.0 90.6 
58.2 36.6 56.8 26.1 
47.7 68.3 85.6 92.8 


TABLE 6 - CENTRIFUGAL AND GRAVITATIONAL LOADS* 


Fy My Fy My F, Mz 

(Ib) (in.-Ib) (Ib) (in.-Ib) (Ib) (in.-Ib) 

Mean, (F, M) 0.036 0.191 0.774 0.0 7.300 -0.510 
First Harmonic Amplitude, (F,M), 0.0 0.0 0.194 0.314 0.244 0.0 
First Harmonic Phase, (Oe, uw), ; 0.0 0.0 -96.0 96.0 -159.0 0.0 

(deg) 

Mean, (F, M) 0.035 0.184 0.744 0.0 6.930 -0.485 

First Harmonic Amplitude, (F,M), 0.0 0.0 0.194 0.314 0.244 0.0 


First Harmonic Phase, (@- m4: 0.0 0.0 -96.0 96.0 - 159.0 0.0 
(deg) ; 


Mean, (F, M) 0.030 0.161 0.642 0.0 5.860 -0.410 

First Harmonic Amplitude, (F,M), 0.0 0.0 0.194 0.314 0.244 0.0 

First Harmonic Phase, (Op mu), ; 0.0 0.0 -96.0 96.0 - 159.0 0.0 
(deg) 

Mean, (F, M) 0.024 u.127 0.493 0.0 4.400 -0.308 

First Harmonic Amplitude, (F,™M), 0.0 0.0 0.194 0.314 0.244 0.0 

First Harmonic Phase, (o-, Mm) 4: 0.0 0.0 -96.0 96.0 -159.0 0.0 
(deg) 

Mean, (F, M) 0.016 0.108 0.434 0.0 3.820 -0.207 

First Harmonic Amplitude, (F,M), 0.0 0.0 0.194 0.314 0.244 0.0 


First Harmonic Phase, (@- m)4. 0.0 0.0 -96.0 96.0 - 159.0 0.0 
(deg) ; 


*The results shown here are for an Aluminum model propeller, Op = 5.44 Ibf-s?/#t? and were 
obtained by fairing the experimental air-spin data using the following restraints: 
(F, M)/n? = constant 
(F,,™M,) =constant 
@,  =constant 


Except where indicated in Table 8, all other centrifugal and gravitational loads presented in this 
report are for a Bronze propeller, Pp = 14.48 Ibf-s2/#t? 
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13.76 


12.65 


10.96 


10.21 


TABLE 7 - SUMMARY OF CIRCUMFERENTIAL VARIATION OF LOADS AT THE 
SELF PROPULSION CONDITION; V = 6.52 KNOTS, 
n = 14.08 REVOLUTIONS PER SECOND 


Maximum Values Minimum Values First Harmonic Values 
(F, M) ax (F,M) vain (ee M), 
a, Tas Puax — Omin ——— ?, 
|F,M| |F,M| (F, M) 
Frey 1.43 124 0.51 308 0.42 116 
My, 1.42 128 0.52 312 0.41 120 
Fv 1.40 120 0.54 292 0.38 105 
My, -1.34* 124 -0.61 300 0.32 294 
Foy -7.88* 252 6.30 68 6.87 54 
zy -1.30* 204 -0.63 64 0.31 41 
Mo.3,, 1.37 132 0.59 312 0.35 124 
Moa 1.36 132 0.60 312 0.34 128 
Au 

F 1.42 124 0.53 308 0.41 116 
My 1.40 128 0.54 Sih2. 0.39 120 
Gy 1.14 132 0.83 304 0.12 123 
M,. -1.26* 136 -0.70 308 0.23 304 
PS 1.04 76 0.96 264 0.03 89 
M, -1.18* 204 -0.77 64 0.19 41 
Mo3 1.45 136 0.48 312 0.43 124 
Mo.4 1.60 140 0.33 312 0.56 126 


*The maximum values shown are the values which have the largest absolute value. For maximum values 
shown with negative sign, the corresponding time average values are negative using the adopted 
convention shown in Figure 1. 
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TABLE 8 - TIME AVERAGE LOADS FOR STEADY-AHEAD OPERATION AT THE 
SELF PROPULSION CONDITION; V = 6.52 KNOTS, 


FL (Ib) 
My (in.- Ib) 
F (Ib) 
M, (in. -Ib) 
F, (Ib) 
M, (in. -Ib) 
M3 (in. - Ib) 
Mo.4 (in. - Ib) 
Mos (in. - Ib) 


For Aluminum blades p, = 5.44 Ibf-s?/ft* 


For Bronze blades pp = 14.48 Ibf - s2/Ft* 


Total with 
Aluminum 
Blades 


3.722 
10.903 
3.456 
-7.765 
7.145 
- 2.608 
6.659 
4.653 
2.790 


Total with 
Bronze 


Blades 
3.781 
11.220 
4.742 
-7.765 
19.279 
- 3.456 
5.563 
3.314 


1.079 


Hydrodynamic 


Loads 
Only 


3.686 
10.712 
2.682 
-7.765 
-0.154 
- 2.098 
7.322 
5.609 
3.820 


Kp = F/(pn2D*) and Ku = M/(pn2D°) are based on the density of water 
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n = 14.08 REVOLUTIONS PER SECOND 


Total with 


Bronze 
Blades 


0.0448 
0.0161 
0.0562 
-0.0112 
0.2253 


Hydrodynamic 
Loads 
Only 


0.0437 
0.0154 
0.0318 
-0.0112 
-0.0018 
-0.0030 
0.0106 
0.0081 


0.0055 
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RENSSELAER/DEPT MATH 
2 
ST JOHNS U 
SWRI 
1 APPLIED MECH REVIEW 1 
1 ABRAMSON 
1 BURNSIDE 1 
2 
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STANFORD U/ASHLEY 
STANFORD RES INST LIB 
SIT DAVIDSON LAB 

1 LIB 

1 BRESLIN 

1 TSAKONAS 

1 VALENTINE 
TEXAS U ARL LIB 
UTAH STATE U/JEPPSON 
U WASHINGTON APL LIB 
WEBB INST 

1 LEWIS 

1 WARD 
WHOI OCEAN ENGR DEPT 
WPI ALDEN HYDR LAB LIB 
ASME/RES COMM INFO 
ASNE 
SNAME 
AERO JET-GENERAL/ BECKWITH 
ALLIS CHALMERS, YORK, PA 
ARCTEC, INC/NELKA 
AVCO LYCOMING 
BAKER MANUFACTURING 
BATH IRON WORKS CORP 

1 HANSEN 

1 FFG PROJECT OFFICE 
BETHLEHEM STEEL NY/DE.LUCE 
BETHLEHEM STEEL SPARROWS 
BIRD-JOHNSON CO 


1 CASE 
1 RIDLEY 


Copies 


BOEING ADV MAR SYS DIV 
BOLT BERANEK AND NEWMAN 
1 BROWN 
1 JACKSON 
BREWER ENGR LAB 
CAMBRIDGE ACOUS/JUNGER 
CALSPAN, INC/RITTER 
EASTERN RES GROUP 
EXXON DES DIV 


FRIEDE & GOLDMAN/MICHEL 


GEN DYN CONVAIR 
1 ASW-MARINE SCIENCES 


GEN DYN ELEC BOAT/ 
BOATWRIGHT 


GIBBS & COX 

1 TECH LIB 

1 OLSON 

1 CAPT. NELSON 
GRUMMAN AEROSPACE/CARL 
HYDRONAUTICS 

1 DUNNE 

1 SCHERER 

1 LIBRARY 
INGALLS SHIPBUILDING 
INST FOR DEFENSE ANAL 


ITEK VIDYA 


LITTLETON R & ENGR CORP/REED 


LITTON INDUSTRIES 
LOCKHEED M&S/WAID 
MARINE VIBRATION ASSOC 


1 BRADSHOW 
1 VASSILOPOULIS 
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Copies 


Copies 


ee 


DOUGLAS AIRCRAFT 

1 HESS 

1 SMITH 
NATIONAL STEEL & SHIPBLDG 
NEWPORT NEWS SHIPBLDG LIB 
NIELSEN ENGR/SPANGLER 
OCEANICS/KAPLAN 
NAR SPACE/UJIHARA 
PROPULSION DYNAMICS, INC 
K.E. SCHOENHERR 
GEORGE G. SHARP 


SPERRY SYS MGMT LIB/ 
SHAPIRO 


SUN SHIPBUILDING 
1 PAVLIK 
1 SINGH 
ROBERT TAGGART 


TETRA TECH PASADENA/ 
CHAPKIS 


TRACOR 


UA HAMILTON STANDARD/ 
CORNELL 


CENTER DISTRIBUTION 


Code 

ILAL Ellsworth 
1102.1 Nakonechny 
15 Cummins 
1509 Pollard 
152 Wermter 
1524 Lin 

1524 Hecker 
1524 Day 

1524 Remmers 
1524 Gordon 


Copies Code 


1 1524 Roddy 

1 1532 Dobay 

1 1532L Denny 

1 154 Morgan 

1 1544 Cumming 

30 1544 Boswell 

10 1544 Jessup 

1 156 Hagen 

all 1556 Santore 

iL 1556 Wisler 

1 1556 Jeffers 

1 LTD Krenzke 

1 1720.6 Rockwell 

iL 19 Sevik 

i 1962 Zaloumis 

iL 1962 Noonan 

il 1962 Antonides 
il PIs) Spargo 

iL 2731 Moken 

1 2814 Czyryca 

30 5214.1 Reports Distribution 
iL 5221 Library (C) 
ilk 5222 Library (A) 
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DTNSRDC ISSUES THREE TYPES OF REPORTS 


1. DTNSRDC REPORTS, A FORMAL SERIES, CONTAIN INFORMATION OF PERMANENT TECH- 
NICAL VALUE. THEY CARRY A CONSECUTIVE NUMERICAL IDENTIFICATION REGARDLESS OF 
THEIR CLASSIFICATION OR THE ORIGINATING DEPARTMENT. 


2. DEPARTMENTAL REPORTS, A SEMIFORMAL SERIES, CONTAIN INFORMATION OF A PRELIM- 
INARY, TEMPORARY, OR PROPRIETARY NATURE OR OF LIMITED INTEREST OR SIGNIFICANCE. 
THEY CARRY A DEPARTMENTAL ALPHANUMERICAL IDENTIFICATION. 


3. TECHNICAL MEMORANDA, AN INFORMAL SERIES, CONTAIN TECHNICAL DOCUMENTATION 
OF LIMITED USE AND INTEREST. THEY ARE PRIMARILY WORKING PAPERS INTENDED FOR IN- 
TERNAL USE. THEY CARRY AN IDENTIFYING NUMBER WHICH INDICATES THEIR TYPE AND THE 
NUMERICAL CODE OF THE ORIGINATING DEPARTMENT. ANY DISTRIBUTION OUTSIDE DTNSRDC 
MUST BE APPROVED BY THE HEAD OF THE ORIGINATING DEPARTMENT ON A CASE-BY-CASE 
BASIS. 


